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43 ■ Abstract 

> : 

^ ■ We show how sum rules for the weak decays of heavy flavor hadrons can be 

derived as the moments of spectral distributions in the small velocity (SV) limit. 
This systematic approach allows us to determine corrections to these sum rules, to 
obtain new sum rules and it provides us with a transparent physical interpretation; 
it also opens a new perspective on the notion of the heavy quark mass. Applying 
these sum rules we derive a lower bound on the deviation of the exclusive form 
factor Fb~>d* from unity at zero recoil; likewise we give a field-theoretical derivation 
of a previously formulated inequality between the expectation value for the kinetic 
energy operator of the heavy quark and for the chromomagnetic operator. We 
analyze how the known results on nonperturbative corrections must be understood 
when one takes into account the normalization point dependence of the low scale 
parameters. Relation between the field-theoretic derivation of the sum rules and the 
quantum-mechanical approach is elucidated. 
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1 Introduction 



The theory of preasymptotic effects in inclusive weak decays of heavy flavor hadrons 
has been developing from the early eighties on — || and is entering now a rather 
mature stage || |19 . In recent papers [15, 16, 17| it was shown in particular 
how the effects due to the motion of the heavy quarks inside the hadrons can be 
incorporated in a systematic way, namely through distribution functions which cru- 
cially depend on the ratio 7 = m q /mQ where rriQ and m q are the masses of the 
initial and final state quarks, respectively. Among other things, it was mentioned 
that the formalism developed in Ref. |15]| automatically ensures the Bjorken sum 
F2D], [2T|, f£J in the small velocity (SV) limit g§ (Sect. 4 of Ref. (I§) 
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In 
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the present work we discuss in more detail the sum rules of Refs 
the so-called optical sum rules derived later by Voloshin p4] (see also [25 
well as other similar sum rules including those considered by Lipkin within a 
quantum-mechanical approach. 

These are the main observations of our paper: 

• We demonstrate that these apparently isolated sum rules represent merely 
different moments of observable spectral distributions. Their physical meaning be- 
comes absolutely transparent within the formulation of the problem suggested in 
Ref. |TJ|. Actually many crucial elements are already included, implicitly and ex- 
plicitly, in Ref. [HJ, so our task is to combine them. This approach will allow us 



to get new sum rules and to obtain corrections to the old ones in a systematic and 
comprehensive way. 

• The sum rules are used to derive a lower bound on the deviation of the ex- 
clusive form factor Fb^d* from unity at zero recoil. The lower bound is essentially 
determined by the average value of the chromomagnetic operator /j,q, familiar from 
previous studies. We also obtain a new field-theoretic derivation of the previously 
formulated inequality /i^ > (Mq, where 
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Its validity has been questioned by some authors under the pretext that the original 
line of reasoning was purely quantum-mechanical. 

• The analysis of the sum rules gives us an opportunity to discuss the notion of 
the heavy quark mass from a point of view complementary to recent investigations 
p8| , p£fl . It is shown that the key theoretical parameter A(/i) is directly related 
to quantities measurable in inclusive semileptonic decays of B mesons in a certain 
kinematical regime. The relation obtained makes absolutely explicit the fact that 
A does not exist as a universal constant, as had been previously believed. Any 
consistent treatment must deal with a fi dependent function, A(/i), where \i is a 
normalization point. 

The assertions formulated above for A(/i) are applicable in full to the expectation 
value of the kinetic energy operator /j, 2 . 
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• The conventional derivation of the sum rules in the SV limit based on QCD 
and the heavy quark expansion is shown to be equivalent to Lipkin's approach |27] 



in all aspects where the ultraviolet domain (i.e. virtual momenta of order of the 
heavy quark masses) is unimportant. In Lipkin's formalism the problem is treated 
quantum-mechanically from the very beginning; the b —>■ c transition is treated as an 
instantaneous perturbation. These two formalisms are complimentary with respect 
to each other. 

Generically the sum rules considered yield the moments of energy of the final 
hadronic state in the b — > c transition as an expansion in different theoretical param- 
eters: the perturbative expansion in a s , the nonperturbative expansion in l/m& )C , 
and, finally, the expansion in q/m c where q is the four-momentum carried away by 
the current (the current which induces the b — > c transition under consideration). 
This latter expansion is nothing else than the SV limit. In the point of the zero 
recoil q = 0. Near zero recoil \q\/m c can be treated as a small parameter and all 
quantities can be consistently expanded in this parameter. Each particular term in 
this expansion generates its own set of sum rules. We will be mostly dealing with 
the sum rules at zero recoil and those emerging in the order (q/m c ) 2 . 

In the focus of our investigation are the nonperturbative corrections proportional 
to powers of l/m^ c . The Wilson operator product expansion (OPE) |30| is the the- 
oretical foundation for calculation of these corrections. The version of OPE used 
in the problem at hand is somewhat peculiar and is supplemented by a systematic 
expansion in the inverse heavy quark masses which is close, both conceptually and 
technically, to what is currently known as Heavy Quark Effective Theory (HQET) 
In many works devoted to the subject it is not realized that HQET can be 
consistently formulated only in the context of OPE, with all attributes needed for 
implementation of Wilson's idea of separation of large and short distance contribu- 
tions. In particular, the theory of the l/m^ c corrections cannot be built without 
explicitly introducing a normalization point /x, the boundary between the large and 
short distance domains. All theoretical parameters appearing in the expansion are 
//-dependent, generally speaking; at the same time, the predictions for measurable 
quantities are //-independent, of course. The issue of the \x dependence is implicit as 
long as one neglects perturbative corrections. Simultaneous analysis of the pertur- 
bative and nonperturbative terms calls for explicit introduction of the normalization 
point. This aspect is discussed in a special section, and we repeatedly emphasize, 
throughout the paper, that all HQET-type expansions become well-defined only 
provided they are formulated in the framework of the operator product expansion. 

The organization of the paper is as follows. In Sect. 2 we outline the specifics 
of the general OPE-based approach applied to heavy quark decays and derive the 
expansion for the heavy flavor hadron masses in QCD. In Sect. 3 basic elements of 
our approach are demonstrated in a simplified model where the heavy quarks are 
deprived of their spins and the external "weak" current considered is scalar. Sections 
4 and 5 are devoted to real QCD; here we derive a set of sum rules describing the 
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inclusive semileptonic decay B — > X c Zz/. The sum rules are considered in detail 
at zero recoil in Sect. 4 and in the general situation with the emphasis on the SV 
kinematics, q -C m c , in Sect. 5. We find a lower bound on 1 — F^ D , at zero recoil 
and present expressions for A(/x) in terms of the differential distributions observable 
in semileptonic B meson decays. In Sect. 6 we establish a relation between field- 
theoretic and quantum mechanical derivations of the sum rules in the SV limit. 
Sect. 7 addresses practical implications of "running" of basic low-energy parameters 
of the heavy quark theory. Sect. 8 briefly summarizes the main results. In Appendix 
we give a quantum-mechanical derivation of a sum rule relating A to observable 
quantities. 

2 General OPE Approach and Mass Formulae 

Inclusive heavy flavor decays are closely related to deep inelastic scattering (DIS). 
While the latter was in the focus of theoretical investigations in the early days of 
QCD, heavy flavor decays received marginal attention. Recently it has been shown 
that many elements of the theory of DIS find their parallels in heavy flavor decays. 
In this paper we discuss one more aspect with an apparent analogy in the theoretical 
treatments, namely the sum rules. 

Let us recall that the standard analysis of DIS |3^] proceeds as follows. One 
starts from the operator product expansion (OPE) for the T product 

%» = i I d 4 x e"^ T {jI(x)UO)} = £ ^^(q) ■ 0« .^(0) (2) 

where j M is the electromagnetic or some other current of interest and 0^(0) are local 
operators. The average of T Mi , over the nucleon state with momentum p presents 
a forward scattering amplitude of the Compton type. This amplitude depends on 
two kinematic variables q 2 and v = qp. For large Euclidean q 2 OPE leads to a 
set of predictions for the coefficients of the Taylor expansion in v of the Compton 
amplitude at v — 0. These predictions are formulated in terms of the expectation 
values of local operators (N\0^\N) over the nucleon state. The coefficients of the 
expansion are related, via dispersion relations, to integrals over the imaginary part 
of the amplitude at hand (moments of the structure functions). 

The strategy, as well as the results obtained, are quite general. At the same 
time, certain moments play a distinguished role due to the fact that they turn out 
to be proportional to operators whose matrix elements between the nucleon state 
are known on general grounds. Relations emerging in this way are called sum rules 
proper and possess particular names. In the case of unpolarized targets we deal with 
the Adler sum rule for neutrino scattering and the Gross-Lewellyn-Smith sum rule 

Conceptually a very similar description can be applied to the weak decays of 
heavy flavor hadrons Hq. One relates the observable quantities to a non-local tran- 
sition operator T, expands the latter into a series of local operators and determines 
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their expectation values between the state Hq. There exist, of course, several tech- 
nical differences relative to the case of DIS: one deals with heavy quark currents, 
uses the heavy quark mass as the expansion parameter instead of the square of the 
momentum transfer, and forms the expectation values for the heavy flavor hadron. 
Moreover, much larger number of these expectation values is known (in the leading 
approximation in l/mg), as compared to light hadrons. It is due to the fact that 
the distribution of the heavy quark inside hadron is trivial to leading order, in con- 
trast with the case of light quarks. As a consequence one can predict, for example, 
absolute decay rates, not only their evolution as the scale changes. 

Now we make a digression of a general nature concerning the heavy quark ex- 
pansion for the hadron masses. The corresponding expressions will appear in many 
instances in the sum rules below. 

For any hadron its mass can be written as a matrix element of the trace of full 
energy-momentum tensor of the theory, in particular, 

Mh q = ^^H Q \9 m \H Q ) y (3) 

where 9^ denotes the energy-momentum tensor; we use the relativistic normaliza- 
tion of the states, 

(H Q \H Q ) = 2M Hq V 

for the hadron at rest. Our goal is to obtain the expansion for Mjj q in the inverse 
powers of rriQ. In particular, we need to construct l/mg expansion for 9^ which 
follows from the expansion of the Lagrangian. 
The original QCD Lagrangian has the form 

£q CD = -^- 2 (G%) 2 + E ™ q )q (4) 

Us q 

where the sum runs over all existing quarks. In the field theory one has to specify 
exactly the normalization point ft where all operators are defined; thus all couplings, 
viz. g s and are functions of fi. The standard Eq. (Q) assumes that the normal- 
ization point is much higher than all masses in the theory, fi ^> rrij. In particular, 
no terms ~ [mj fi) n are kept. 

Constructing the effective theory designated to describe the low energy properties 
of heavy flavor hadrons we need to have the normalization point fi below mass of the 
heavy quark tuq. This changes the generic form of the Lagrangian and a series of 
operators of higher dimension can appear. Their coefficients contain inverse powers 
of the heavy quark mass tuq. It is important that all these effective operators are 
bound to be Lorentz scalars. 

In what follows we are not interested in terms of order higher than l/nig. To 
this accuracy the general form of the Lagrangian is 

C = ~^(G%) 2 + Q(iP- m Q )Q + -£q(ip- mq ) q + -^-Qia^Q + 

^9s q 4m Q 
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+ —T m qQ i(J ^ G ^Q + J2 -^-QTQqTq + J2 -^q^qq'Tq' 




q m Q T,q m Q 

+ Zw2 Tr G^ v G vp G Pil + O ( ^g- ) (5) 

m Q 

where G a v is the gluon field strength tensor and G pv = t a G a v is its matrix rep- 
resentation (Trt a t 6 = 5 a b/2); below we will often use the short-hand notation 
iaG = ia pu G pu = ij^jvG^. The sum over light quark flavors is shown explic- 
itly as well as the sum over possible structures F of the four fermion operators. 
Again, all masses, couplings and coefficients of higher dimension operators depend 
on the renormalization point. For example, the coefficient cq has the form 

_ 3 

^feS)) 1 - 1 ^ • 6 = u -1"' • < 6) 

Expansion in the parameter I/tuq leads to the fact that the logarithmic evolution 
of operators generally mixes them with the ones of higher dimension, in contrast 
with the expansion in m q for light quarks. The above relation is the effect of such 
mixing of the operator QQ with the chromomagnetic one, QiaGQ. 

Evolving the parameters of the effective Lagrangian to the normalization point 
ix one can use the values of c, d, f, h obtained at \l ~ mq as the initial conditions. 
In fact all operators with d > 5 appear at this scale only as a result of quantum 
corrections and contain explicit powers of a s (mg)/7r. For the sake of simplicity we 
will generally neglect them in our subsequent analysis, although it is easy to keep 
track of them if necessary. This neglection does not mean that we disregard all 
effects of the order of \jw?n\ they appear from the leading "tree level" operators as 
well and will be accounted for consistently. Note that neglecting higher order terms 
in Eq. (|5|) gives us the standard QCD equations of motion for both the gluon and 
quark fields. If we were interested in quantum corrections at the order of 1/mq or 
1/tbq we would need to incorporate additional terms in the equation of motions or, 
alternatively, consistently consider these higher order terms in the Lagrangian as 
perturbation and expand in them to necessary order. 

In the theory with Lagrangian Eq. (|4]) the trace of the energy-momentum tensor 
has the form 

W = C + (^TT. ~ m «d^) C ( 7 ) 

m Q QQ , (8) 




(9) 
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where 7 m is the anomalous dimension of light quark mass, //dm ? (/i)/d/i = —7. 
and P(a s ) is the Gell-Mann-Low function, /ida s (/z)/d/i = (3{a s ). Note that the 
scale dependence of mass enters 6^ in the same way for light and heavy quarks; 
the explicit form of this dependence is, of course, different. The modification of 9^ 
for the case when higher dimension operators are included is straightforward using 
the general expression Eq. (^). In Eq. T> has the meaning of the part associated 
with light degrees of freedom. 

In the case we consider, thus, the hadron mass can be expressed in terms of two 
expectation values, 

Mh * = ^—(H Q \m Q QQ\H Q ) + -±-(H Q \V\H Q ) . (10) 



The advantage of this formula over a more familiar representation [31] where the 
mass is calculated through the expectation value of the Hamiltonian is rather trans- 
parent - mass is a Lorentz scalar quantity, and in Eq. (|10j) it is given as the 
expectation value of the Lorentz scalar operators. Of course, if one wants to develop 
1/tuq expansion for M^ Q one eventually arrives at the standard expansion involving 
the same nonrelativistic operators. All relevant terms in the effective Hamiltonian 
appear to be directly related to the 1/tuq expansion of the operator QQ. Let us 
elucidate this point in more detail dealing with the two operators in Eq. (|T0|) in 
turn. 

To calculate the matrix element of QQ we, following Refs. || [7J, start from the 
heavy quark current. In the rest frame of Hq we have: 

' (H Q \Q l0 Q\H Q ) = 1 . (11) 



2M Hq 



Using the decomposition 

iD„ = m QVfl + tt m , VlI = (phq)^ I M Hq (12) 
and the equation of motion i p Q = tuq Q, which imply 

we arrive at the following identity 

QQ = QloQ + 2q( 1 -^-) 2 Q = Q l0 Q - 2Q-^---^Q = 

V 2 J 2m Q 2m Q 

- f 

— QloQ + Q 2~Q + total derivative ; (14) 



7 



in the first relation above the operators 7r M act on the Q field. We are considering 
the forward matrix elements with zero momentum transfer and thus can drop all 
terms with total derivatives. Eq. flH| ) supplemented with equations of motion ([13]) 
generates then the complete 1/mq expansion for the scalar density: 



i 

2m A n ~" N " ' 2 ' ' ' * '" * 2m 



QQ = QloQ + 7^tQ(k 2 + -ctG)Q = Q 7o Q - 7r ^-Q (7?a) 2 Q 



Q 

Q(-{DE) + 2a-Exn)Q + O [ ^j] ■ (15) 




Here Ei = Gio is chromoelectric field and its covariant derivative is defined asQ 
DjEk = — i [ttj , Ek] ; we have omitted the term Q([iTk, [^o,^i\] — [fti, Wo^k]])Q using 
the Jacobi identity. Moreover, DE = g 2 t a J^ by virtue of the QCD equation of 
motion, where J® = J2qQ1^t a( l is the color quark current; therefore the first of the 
1/rriQ terms can be rewritten as the local four fermion interaction [15|. 
We see that the itiqQQ part of Q m in Eq. flTD| ) generates 



mQ -2M^ Q {HQlQ 



-A (-(^ + 2^.^X7? 



Q\h q ) + o\- 



2m,Q 4m^ 

"(16) 

Eq. (|T6|) shows the explicit dependence of 6* MM on the heavy quark mass. However, 
it is not the only source of ttiq dependence. Indeed, the hadronic state Hq also 
depends on ijiq. To account for this dependence one needs to introduce the basis 
of asymptotic (mq = oo) states, \HQ) mQ=00 and develop the perturbation theory in 
explicit 1/rriQ terms in the Largangian. Let us emphasize that the explicit 1/toq 
terms above come from short distances of order I/ttiq whereas effects due to mq 
dependence of states Hq are large distance (~ Aqq D ) ones. 

We are interested in terms ~ Aqcd, Aq CD /mg and AQ CD /mg in mass. Therefore 
we must account for 1/mq effects in the matrix elements of Q(aif) 2 Q as well as terms 
mq 1 and mq 2 in the expectation value of T>. 

The term mq in Eq. (|T6|) is the only one in Mh q which is linear in mass. 
To zeroth order in mq the mqQQ part gives no contribution and the only mq- 
independent term comes from the expectation value of T>. To leading order in mq 
this expectation value presents what is usually called A, 



1 , „ , x , „ , dmr 



(Hq\ V \Hn) = ——(Hq\ ^^-G 2 \Hq) + 



+ ,^-^HQ\m q {l + lrn )qq\HQ) - // ^ = A + O^ 1 ) . (17) 

In a sense, the first term is an analog of the gluon condensate in the QCD vacuum 
P3fl . Instead of the QCD vacuum we deal now, however, with the ground state in 
the sector with the heavy quark charge equal to unity. 

1 Note that in our notation D = —d/dx — iA, therefore (DE) = — divE in the Abelian case. 
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The analogy continues further and allows us to derive low energy theorems very 
similar to those taking place |3|] for the vacuum correlation functions involving G 2 



and other operators. In the present context the low energy theorems take the form 
-i J d A x -J—(H Q \T {v(0),Q(x)OQ(x)}\H Q ) = 



d-^-(H Q \QOQ\H Q ) (l + 0(m c 



;is) 



where the operator Q(x)OQ(x) is bilinear in Q and Q and O contains derivatives 
and light fields A^, q; the factor d is equal to the dimension of O. Derivation of 
these low energy relations merely parallels that given in Ref. | ]34| | . 

Armed with these theorems we will be able to calculate subleading terms in Eq. 
(|17|). But first we need the expansion of the QCD Lagrangian to the necessary order 
in 1/itlq. We obtain in the standard way 



Q(iW-m Q )Q 




1 



710 



2m Q 



7TO- 



8mj 



(-(DE) 



2(7 • E XT? 



0"7T 



5m, 



where 



2 \ rrin 



1 




(19) 
(20) 



is nonrelativistic two-component spinor field. This substitution is nothing but the 
Foldy-Wouthuysen transformation which is necessary to keep term linear in ttq in 
its canonical form to ensure the mass independence of the field (p. If it is not 
accomplished there would be implicit dependence of the heavy quark fields on mass, 
and the low energy theorems will be modified. In the last equation (ITT 



n 



i 



i 



2m f 



0~7T 



-51717 



-(DE) + 2a ■ E 



X 7T 



(«77T) 2 = 7T 2 + OB (21) 



is the nonrelativistic Hamiltonian through second order in I/tuq. It coincides with 
the text-book expression |35| in the Abelian case. Notice that 1/m.g and 1/m.q 
terms in Eq. ( P~5| ) are just dri/dniQ. 

In fact it includes the anomalous dimension of QOQ which appears when radiative corrections 
are incorporated. 
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Now we can calculate power in I/tuq corrections to matrix elements. Indeed, 
the next-to-leading term in {2M Hq )~ 1 (Hq \ V \Hq) is given by 



m Q' 



2M 



H, 



{Hq\V\H q ) 



2 x 



1 / m Q 

f L 



% Id x 



2M, 



H, 



(H Q \T{V(0),Q(x) 



I (JIT 



-(h q \q [ -^Lq\h q ) 



Combining now Eqs. (^) and (|T5|) we arrive at a standard formula pi 
hadron mass, 



Mh q =m Q + A + 



m Q \2M HQ S QIV 2 Vl Ql 



0(m- c 



(22) 
7[ for the 

• (23) 



?TtQ=00 



Using the same technology it is not difficult to develop the expansion one step 
further to include the 0(rriQ 2 ) term. To this end we must find and take into ac- 
count the 1/rriQ contribution in (Hq\Q(<j7t) 2 Q\Hq) in Eq. flT5|) as well as the I/tuq 
contribution in (Hq\G 2 \Hq) . For the first matrix element we have 



/ — 



Hq 



■(H Q \Q(at) 2 Q\H Q ) 



1/m.Q 



i Id x 



AM. 



-(H Q \T{Q(x)(a^ 2 Q(x),Q(0)(^) 2 Q(0)}\H Q y r 



m Q -- 



where p 3 is a positive parameter of the order of Aq CD measuring the correlation 
function above. The prime in Eq. (|24]) indicates that the diagonal transitions 
have to be removed from the correlation function (analogously to elimination of the 
disconnected parts in the vacuum correlators). 

By the same token using low energy theorems, the l/rrig piece in (2M# Q ) -1 - 
(Hq \ T> \ Hq) can be written in terms of a similar correlation function 



rri 



Q 



2M 



Hq 



(Hq\V\H q ) 



l/m 2 



2M 



-I] 



Hq 



1 = 



d*x (Hq\T V(0), -Q(x) {-(DE) + 2a-Exn) Q(x) \H Q ) - 



8 2M 



Hq 



d A xd\ (H Q \T{v(0),Q(x)(an) 2 Q(x),Q(y)(a7c) 2 Q(y)}\H Q ) 



8 2M Hq 



(H q \ Q(x) (—(DE) + 2d ■ E x 7?) Q(x)\H Q ) 



(25) 
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where the proper analog of Eq. fll8|) is used. 

Now we have at our disposal everything needed to write down the mass formula 
including 1/toq corrections. Combining Eqs.(|Tt3|), (^|) and ( [25] ) we obtain 



[say 



M HQ = mQ + A + -L {-J_ (ff0 | Q ^ Q |i/ 0> > + 



rriQ=oo 
o3 



(26) 

It is clear that the rriQ dependence obtained has the usual quantum mechanical 
interpretation: the third and fourth terms in Eq. fl26| ) present the expectation values 
of the Hamiltonian TC, Eq. fl2~T|) , over the asymptotic state \Hq) mQ=00 , and the term 
— p 3 /4mg presents the second order iteration of the Hamiltonian. 

For further usage we introduce here the following notations for the expectation 
values of the two terms in 1/tuq Hamiltonian, which are the Darwin and the con- 
vection current (spin-orbital) interactions, respectively: 

As = ^{B\bS^Exnb\B) = (B'\b 9 ■ E x 9 b\B') . (27) 

p 3 D is directly related to the third moment of the heavy quark distribution function 
and has been estimated in Ref. [15| (see also [36|). On the other hand, the value of 
p\ s is expected to be suppressed for L = states like B and B*. 

A similar decomposition can be obviously made for the non-local correlators 
whose expectation value has been generically denoted by p 3 ; the four unknown 
parameters appearing here can be related to the 1/tuq part of the expectation 
values of Qti 2 Q and QaBQ in pseudoscalar and vector ground states: 

fL = if d'x ^-(B\T{bn 2 b(x), b^b(0)}\BY 
p " G = t I * x ^F< B I T ^ 26 W' bffBb(0)}\BY 



TjplMfc* + -PaWh ~ S a 5 jk ) =i [d A x —-(BlTibaiB^x), ba^b^B)' . 

; (28) 

Then one has for the parameters p 3 in B and B*, respectively 



( p 3) B = P l 7T + P l G + P l + p \ , ( p 3 )^ = P l n - Ip 3 G + p 3 - Ip 3 . (29) 
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The l/rriQ corrections to the expectation values of tt 2 and aB are 



2M R " ,d ™ 2m h 



1 l->2\ 2 2P7T7r 3P7rG i — 2^ 



2M B * n ' r ^ 2m fe 

- A*G o h U ( m b 



2Mb 2m b 

.3 i «„3 o„3 



1 (aB) B * = \& - ^±gg 2 ^ + 0(m 6 -) • (30) 



2Mb* "i Gm b 



Very similar parameters have been introduced in Ref. |36| where the mass formulae 
were discussed in the standard HQET approach 

In the subsequent section we will consider in detail the toy model where heavy 
quarks are spinless and light quark masses vanish [IS] . The absence of spin simplifies 
the analysis above. Let us briefly review the changes. In this toy model the trace 
of energy-momentum tensor is given by 

»~ = a "^ + S G *~ G *~" 2 ^N 0i (31) 

the heavy quark field is still denoted by Q but its dimension is m now, not m 3//2 . 
Moreover, this expression, as well as the most of other ones can be obtained from 
the spinor case by substitution 



Qspinor * y Qscalar • 

Spin matrices are now absent, i.e. (<?7?) 2 p inor — * (7?) 2 C aiar- I n particular, the equation 
of motion and the expansion of the expectation value for QQ take the form: 

7T 2 

n Q = -- — Q 

-±—(H Q \2m Q QQ\H Q ) = 1 - -±—(H Q \Q— Q\H Q ) + 0(m^) . (32) 
Hq zm h q ™q 

Notice that in the scalar case there is no explicit 1/rriQ terms; it corresponds to the 
absence of 1/?tjq terms in the nonrelativistic Hamiltonian for scalar particles: 



3 We do not agree with some numerical coefficients of Ref. |3(|. In particular, Eqs. (50-55) of 
that paper are not quite consistent. Also, our result for the contribution of the Darwin term in 
mass is positive in the factorization approximation and smaller by a factor of two. 
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Low energy theorems have exactly the same form as before with 



As a result the expansion for the hadron mass in the scalar case takes the form 

m h q =m Q + A + —J— {(H Q \Qn 2 Q\H Q )} - ? + 0(m Q 3 ) (33) 



Q 



with 



iriQ- 



p 3 = 2m 2 Q i J d 4 x^-(H Q \T{Q(x)7r 2 Q(x),Q(0)7r 2 Q(0)}\H Q y (34) 



and 



(Hq\Qti 2 Q\Hq) = ¥*L {(H Q \Qn 2 Q\H Q )\ - + 0{m Q 2 ) . (35) 

rriQ ' m Q =oo rriQ ■* 

Let us return to ordinary QCD and briefly discuss the normalization point de- 
pendence of A. In this context we can neglect terms of order 1/m.q and higher; 
Eq. fllTD gives the following definition of A(/i): 

10.) - {M„ Q - m Q M} mQ _ = ^i-(H Q |^G 2 +X>,U + 7„)?5|//«) ~ 

- ■ (36) 

d/i 

The last term is specific for the field theory and is absent in the naive quantum 
mechanical approach. 

In the perturbation theory the /i-dependence of uiq appears in the order a s , 

— = -Cm h 0{a s ) . (37) 

The matrix element in the right hand side of Eq. fl36| ) constituting the first part of 
A, does not undergo perturbative renormalization f] in order a s . Nevertheless it is 
//-dependent in the theory with heavy quarks; its scale dependence appears in terms 
~ a 2 . It might be tempting to drop the last term in Eq. (p6|) and consider 

5 = 2M^ {h ^tB, g2 + + = Mh « " mM + (38) 



4 The term with the anomalous dimension of the light quark mass per se is not renorm-invariant 
in the order a 2 ., but its /^-dependence is compensated by the mixing with the G 2 term. 
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as an option for A. However, it would not help defining a "purely nonperturbative" 
A; in particular, an attempt to put fi —>■ would lead to the same infrared renor- 



malon as in the pole mass |f28|, |29f . The transparent physical interpretation of this 
phenomenon will be discussed below in Sects. 3.4 and 5.1. Though it is worth noting 
that the combination mg(yu) — /idmg/d/i in Eq. Q3"8| ) has the meaning of the one 
loop pole mass of the heavy quark and therefore it naturally enters the "practical 
OPE" calculations performed at the one loop level (for the detailed discussion see 
Sect. 7). 

The 1 / tuq terms in mass Mh q obtained in this section will be used in derivation 
of the second sum rules at zero recoil. Eq. (|36| ) might be interesting by itself, 
though, as an alternative - compared to the standard HQET analysis - definition 
of one of its key parameters, A. 

3 The Bjorken, 'Optical' and Other Sum Rules — 
Toy Model 

To introduce the reader to the range of questions to be considered below in the most 
straightforward way we first "peel off" all inessential technicalities, like the quark 
spins, and resort to a simplified model which has been previously discussed in Ref. 



1 15 1 - It will be a rather simple exercise to return afterwards to standard QCD which 



will be done in Sect. 4. 

3.1 Description of the model 

We consider a toy example where all quarks are spinless; two, denoted by Q and q, 
couple to a massless real scalar field 0: 

= hQ<Pq + h.c. , (39) 

where h is the coupling constant and Q = . The masses of the quarks Q and q 
are both large. Later on we will analyze the SV limit where 

Aqcd < rn Q - m q < m Q (40) 

has to hold, but for the time being the ratio 7 = m q /mQ can be arbitrary provided 
that mQ—m q ^> Aq CD []. The field carries color charge zero; the reaction Q — > q+<p 
is thus a toy model for the radiative decays of the type B — > X s ^. 

The total width for the free quark decay Q — > q + (f) is given by the following 
expression: 

r(Q - #) = ^ = r (41) 



5 To reach the SV limit in the B — > X c ev transitions there is no need to assume that mi, — m c <C 
ra c . The small velocity regime for the c quark can be ensured by adjusting q 2 appropriately. 
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where 



2 2 

Eo = q ~- 42 

2mn 



As explained in Refs. U ]l5j the theory of preasymptotic effects in inclusive 



decays is based on introducing the transition operator, 

f = ijd i x e- iqx T{Q(x)q(x) , <?(0)Q(0)}. (43) 

Then the energy spectrum of the <fi particle in the inclusive decay is obtained from 
T in the following way: 



dT h 2 E 



-lm(H Q \f\H Q ) (44) 



dE Att 2 M Hq 

where Hq denotes a hadron built from the heavy quark Q and the light cloud 
(including the light antiquark). If not stated otherwise, Hq will denote the ground 
state in a given channel. Moreover, one can (and must) apply the Wilson operator 
product expansion ]30| (OPE) to express the non-local operator T through an infinite 



series of local operators with calculable coefficients. 

3.2 OPE and predictions for observable quantities 

In the Born approximation the transition operator has the form (Fig. 1) 



T = - / d*x(x\Q— — -y -Q\0) , (45) 

J (P - q + 7TJ - m A q 

iDf, = (P ) M + ?V = m Q v [i + ?V 

which is particularly suitable for constructing the OPE by expanding Eq. ( fiEf) in 
powers of 7r (see Ref. [|15] where all notations have been introduced). The operators 
appearing as a result of this expansion are ordered according to their dimension. 
The leading operator QQ has dimension 2 (let us recall that the scalar Q field has 
dimension 1 in contrast to the real quark fields of dimension 3/2, which leads in 
particular to different normalization factors; we still use relativistic normalization 
in the bulk of the paper, except Sect. 6). Its expectation value for the state Hq 
reduces to unity to leading order in I/tuq] this contribution gives rise to the parton 
result ([|1]). Beyond the leading approximation, according to Eq. (|5^), it takes the 
form 

-L-(h q \qq\h q ) = J_ (l - _i 

2M Hq 2m Q \ 2m Q 



-(Hq\QQ\H q ) = — ( 1 - —(H Q \Q^Q\H Q ) + ...); (46) 



thus, one gets a correction of order 1/mq- 

Corrections of the same order come from higher-dimensional operators in the 
expansion of the transition operator. There are no relevant operators of the next- 
to-leading dimension three (more exactly, as first noted in Ref. Q in the framework 
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of HQET, they vanish because of the equations of motion []). The only operator 
of dimension four in our toy model has the form Qt7 2 Q, and it generates a 1/ttiq 
correction after taking the matrix element over Hq. 
After some simple algebra one finds 

vr W \2m Q 12m Q ) 1 0> 

where operators of higher dimension have been ignored, and we have used a short- 
hand notation for the expectation value over Hq: (...) = (Hq\...\Hq) . 

The expansion of ImT into local operators generates more and more singular 
terms at the point where the spectrum would be concentrated in the free quark 
approximation. The physical spectrum on the other hand is a smooth function of 
E. In principle, one could derive a smooth spectrum by summing up an infinite set 



of operators to all orders (for more detail see Refs. [|T^, [16], jlTI). There is no need 
to carry out this summation here, however, since we are interested only in integral 
characteristics (we will discuss certain sum rules), and as far as they are concerned, 
the expansion in Eq. pT| ) is perfectly legitimate. 

At first, we calculate the total width by substituting eq. ([47]) into eq. ([44D and 
integrating over E. As a matter of fact, the result has already been given in Ref. 

HI, 



dE JE M 1 -^ ") (48) 

where the integration runs from to the physical boundary Eq Hvs , expressed in the 
hadron masses 

A/f 2 _ /If 2 

E phys = __Q q_ M M , 4g) 

and the same convention for M q . We use here the standard notation for the expec- 
tation value of 7? 2 , 

ti = 7^(H Q \2m Q Qn 2 Q\H Q ) . (50) 
zm h q 

The expression in the brackets in Eq. (|48| ) is nothing else than the (corrected) 
matrix element of the operator QQ; the only possible effective operator Qtt 2 Q of 

6 This statement is sometimes erroneously interpreted as a proof of the absence of a term linear 
in 1/toq in the total width, see below. As a matter of fact the authors of Ref. j| believed that 
the linear term may appear in the total width from the overall normalization of {HQ\h v h v \Hq) , 
as explicitly stated, e.g., on page 404 of Ref. Q). Also, "matching" to QCD, the necessary step in 
the HQET approach, was not considered. In fact, the question of absence or presence of the 1/rriQ 
correction cannot be solved in HQET per se since the problem requires a full-QCD analysis of the 
total decay rate to perform such a matching, see the corresponding discussion in Ref. |2§] . 
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dimension 4 is not a Lorentz scalar and thus cannot appear in T |L5| . The meaning 
of this term in (QQ) is quite transparent: it reflects time dilation for the moving 
quark, and the coefficient (—1/2) could, therefore, have been guessed from the very- 
beginning. The absence of a correction of order I/tuq in the total width in the toy 
example under consideration is a manifestation of the general theorem established 
in Refs. || |7| and discussed in more detail recently in Ref. |28| . 



As we will see shortly, Eq. fl48|) treated in the SV limit is equivalent to two 
results simultaneously: that of Ref. [^J and the Bjorken sum rule |20|, |21|, with 



appropriate corrections due to terms which were not considered in Refs. 21, 23 



All sum rules analogous to that of eq. ( f48|) below will generically be referred to as 
the first sum rule. 

Eq. (^) is the first example of the sum rules we will be dealing with throughout 
the paper. Its derivation (as well as that of all similar sum rules) intuitively is 
perfectly clear - the integrated spectrum of the decay obtained at the quark-gluon 
level using the OPE is equated to the integrated physical spectrum, saturated by the 
genuine hadronic states. We will elaborate on the justification for this procedure in 
Sect. 4.1 where, among other things, we discuss the accuracy one may expect from 
relations of this type; the general idea lying behind all such relations is wide-spread 
in QCD. For the moment we just adopt the heuristic attitude outlined above without 
submerging into further, less pragmatic, questions. 

Next, we calculate the average energy of the <fi particle. The corresponding sum 
rule in the SV limit is just a version of Voloshin's optical sum rule! 

To be more exact, let us define moments I n 

f E phya i rfr 

'■ = jt dE{E '' h "- E) "r iE < 51 > 

and consider the first moment I\. Notice that in the SV limit Eq ys — E reduces to 
the excitation energy of the final hadron produced in the decay. The factor Eq Hvs — E 
in the integrand eliminates the "elastic" peak, so that the integral is saturated only 
by the inelastic contributions. Let us note in passing (we will discuss this point later 
in more detail) that the optical sum rule [24| in its original formulation is actually 
divergent; a cutoff must be introduced by hand. In our formulation, of course, 
there is no place for divergences; the expectation value of Eq Kvs — E is defined via 
a convolution with the physical spectrum, and since the <fi energy in the decay is 
finite, the expectation value of Eq ys — E, as well as higher moments, are certainly 
finite. The decay kinematics provides us with a natural cutoff at the scale of the 
energy release. 

Explicit calculation of I\ using Eq. ([47]) yields 

,.2 pphys 

Ix = A-^V (52) 

where A is defined as follows 

A = El hys - Eq, (53) 
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see Eqs. (03) and (f49|). The sum rule (|5^) contains Voloshin's result, again with 
corrections left aside previously [p4 , p5f| . The sum rules analogous to that of Eq. 



(p>2|) below will be generically referred to as the second sum rule. 

It is quite evident that the series of such sum rules can readily be continued 
further. For the next moment, for instance, we get 

r E phys 1 AT a 2 P 2 

I 2 = l' dE{ ^-Efl£ = ^ + ^. (54) 

Relations of this type will be referred to as the third sum rule. Analyzing this sum 
rule in the SV limit one obtains, in principle, additional information, not included in 
the results of Refs. [H], It is worth emphasizing that in Eqs. 



and (|54]) we have collected all terms through order Aq CD , whereas those of order 
Aq CD are systematically omitted. Predictions for higher moments would require 
calculating terms C(Aq CD ) and higher. 

3.3 The SV limit and the sum rules 

We proceed now to discussing the sum rules (|48|), (|52| ) and ( |54]) in the SV limit, i.e. 
in the limit of small velocity v = —q/E where — q and E are momentum and energy 
of the final hadronic system. The assumption Eq. fl4"U| ) means that \v\ <C 1. The 
availability of the extra expansion parameter, \v\ <C 1, makes the SV limit a very 
interesting theoretical laboratory. We have already mentioned that in this limit the 
Bjorken sum rule relates the distortion of the "elastic" peak to the integral over the 
inelastic contributions. Here we will elaborate on this issue. First, we briefly recall 
what is already known and then present new results. 
For technical reasons we will assume that 



Aqcd < fn Q - m q < ym<2A QCD . (55) 

The inequality on the right hand side is not essential and can easily be lifted; it helps, 
however, to make all formulae more concise, and we will accept it for a while, thus 



replacing Eq. (^) by the stronger condition of Eq. (f)5|). We now supplement the 
expansion in Aqcd carried out above by an expansion in v. The natural hierarchy 
of parameters in the domain (|55D is as follows 



v 2 , |w|(A QCD /m Q ), A 2 QCD /m 2 Q . 

Terms of order v 4 -C AQ CD /mg will be omitted. 

Expanding the hadron masses in terms of the heavy quark masses tuq and m q 
one finds for scalar quarks 

2 2 

M Q = m Q + A + £^- + 0(l/m 2 Q ), M q = m q + A + ^- + 0{l/m 2 q ) . (56) 
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Eq. fl5"6|) implies, in turn, that 



A = E?" - E = -A (a + ^Mj - „„ J- + ... (57) 

where for convenience we have introduced an auxiliary parameter t> , 

v = (Mq - M q )/M Q ~ (58) 

This parameter approximately coincides in the SV limit with the velocity of the 
heavy hadron produced in the transition Q — > qcf). Indeed, if the mass of the 
produced excited hadron is M q + e then its velocity is 

l -l e v ° e 

I s " = "» " w Q ~ w Q < 59 > 

plus terms of the higher order in e and/or Vq. 

With all these definitions the sum rules (|48|), (|52f) and take the form 



7 = 1 -J| + ... , (60) 

/l = 5^- t '»^ + -' (61 » 

h = \vlvl + ... (62) 

where the ellipses denote the systematically omitted terms of order Aq CD , as well as 
(9(t>QAQ CD ) terms for I x and 0(fgAQ CD ) ones for I 2 . The first term on the right hand 
side of Eq. ( |6lD corresponds to Voloshin's relation, the second term is a correction 
whose physical meaning will soon become clear. 

To interpret the sum rules derived it will be instructive to consider the transition 
operator off the physical cut. Expressing Eq. ([IT]) in terms of E — Eq Hvs we then get 




M^(H Q \f\H Q ) = — \( 1 - J^- ) + 



3mlJ E-El' 

^ Qm 2 Q ) (E-El hy y 3m 2 Q (E - E% hys ) 3 "J ' 1 ' 

To zeroth order in v or Aqcd there remains only the first term on the right 
hand side of Eq. d63|) . The inclusive decay rate is then totally saturated by a 
single "elastic" channel, the production of the ground state meson containing q. 



This is the perfect inclusive-exclusive duality noted in Ref. [23]. The peak in the 



spectrum obtained in the quark transition Q — > (f>q survives hadronization in this 
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approximation; at the hadronic level Hq — > we still have the same peak at 
the same energy. Eqs. (|61|) and (|62D are, of course, trivially satisfied in this case 
because in the absence of the inelastic contributions both sum rules yield vanishing 
numbers. 

Furthermore, the terms with Aqcd come with v. To order A there is only one such 
term, appearing in eq. (|61|) for /j. This term shows that the inelastic production 
must already be present at this level. It corresponds to the production of a meson 
H* with excitation energy ~ Aq CD and a residue oc v^. Then eq. ( j60| ) implies that 
the height of the elastic peak is reduced by Vq. A rough model exemplifying this 
picture can be obtained from Eq. In this approximation, one can rewrite it as 

follows, omitting all numerical factors: 



M oW H Q} = ^ |(1 - ^—^ + < {E _ Et ,, + - K) + -..) • (M) 

The first term is the elastic peak while the second is an inelastic contribution. Of 
course, Eq. (|64| ) is an illustration and not a unique solution. 

Technically, the term v 2 A in Eq. (|6ll) arises because the "elastic" peak of the 
quark transition situated at E — Eq is slightly shifted when we pass to the hadronic 
transition; the genuine hadronic elastic peak is situated at E — Eq Hvs , to the right 
of the quark peak (see Fig. 2). 

Including 0(A ) contributions we get correction terms in Eqs. (|60| ) and (|61|), 
and Eq. (^) for I 2 becomes non-trivial. The new term in Eq. (^) has a simple 
meaning. In the toy example at hand the excited mesons produced in the decay 
at order v 2 are spin-1 mesons, with the vertex proportional to (ve) where e is the 
polarization vector and v is the velocity of the given meson. This velocity, differs, 
however, from vq by terms of order Aqcd/^q, see Eq. fl59l). This rather trivial 
shift in velocity nicely explains all the sum rules above. Indeed, let us take into 
account that the physical v 2 is reduced by the amount ~ v (A/Mq). The height of 
the inelastic contribution is proportional to the square of the physical velocity, while 
the size of the inelastic domain is ~ A. Hence, with our accuracy the right-hand 
side of Eq. ([H]) is expected to be ~ v^A — t>o(A /tuq) , while the right-hand side of 
Eq. (^) is expected to be ~ A t> 2 , in full accord with what we actually have. A 
qualitative model of saturation now takes the form 




(65) 



where 



X 2 /4 A 

A - -j- ~ Aqcd 
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3.4 Perturbative gluon corrections 

So far we have assigned the gluon field to play the role of a soft medium to incorpo- 
rate the effects of long distance dynamics and have completely ignored perturbative 
gluon corrections. Yet those have to be included; among other things the emission 
of hard gluons generates the spectral density outside the end-point region which is 
very relevant for our analysis. 

In calculating radiative gluon correction we can disregard, in the leading ap- 
proximation, nonperturbative effects, like the difference between wiq and Mq or the 
'Fermi' motion of the initial quark. Thus we deal with the decay of the free quark 
Q at rest into q + <j)+ gluon. The virtual gluon contribution merely renormalizes the 
constant h in the analysis presented above. The analogous renormalization for the 
spinor case has been calculated in Ref. [23]. 



The effects from real gluon emission are most simply calculated in the Coulomb 
gauge, where only the graph shown in Fig. 3 contributes. A straightforward com- 
putation yields to leading order in v 2 

dT^ _ 8a s E 3 1 

dE 9tt E m 2 Q Eq — E 1 ' 

It is well-known that the logarithmic singularity in the integral over E for T tot is 
canceled by a contribution of soft virtual gluons to the renormalization of h. For 
the second sum rule we are going to discuss here, this infrared range is not singular. 

Gluon emission obviously contributes to the spectrum in its entire domain < 
E < Eq. In this order a s does not run, of course. To demonstrate its scale de- 
pendence one has to carry out a two-loop calculation; it is quite evident, however, 
that it is a s (E — Eq) that enters. Therefore, strictly speaking, one cannot apply Eq. 
(H|) too close to Eq. Even leaving aside the blowing up of a s (E — Eq), there exists 
another reason not to use Eq. (|66| ) in the vicinity of E : HE is close to E , the 
emitted gluon is soft; such gluons are to be treated as belonging to the soft gluon 
medium in order to avoid double counting. The separation between soft and hard 
gluons is achieved by explicitly introducing a normalization point fi. The value of 
\x should be large enough to justify a small value for a s {p). On the other hand 
we would like to choose \i as small as possible. The possible choice is to have /i 
proportional to Aq CD , but with a constant of proportionality that is much larger 
than unity 

li = CAqcd, C > 1. 

Then we draw a line: to the left of Eq — \i the gluon is considered to be hard, to 
the right soft. Of course, the consistent introduction of the infrared renormalization 
point n requires that the purely perturbative corrections to the weak vertex, even 
if they happen to be infrared convergent, have to be calculated using this explicit 
cutoff as well (see, e.g. the discussion in Ref. f2~8fl ) - which is almost never done in 
practice. The corresponding modifications will be discussed below in Sect. 7. 
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Let us discuss now the sum rule corresponding to the first moment of E Q hys — E, 
i.e. an analog of Eq. ( |BT]) with radiative corrections now included. Since our main 
purpose in this section is methodical, we will limit ourselves to the first order in A 
and the second order in v. A qualitative sketch of how dT /dE looks like is presented 
in Fig. 4. Then the prediction for 1% can obviously be rewritten as 

h = l° { E^-E)-dE = 



V 1 2 

where by definition 



- ., .^o-M I6a s E 3 IT " 



(67) 



A(u) = / " ^^-(El hys - E) dE. (68) 

Without the radiative tail the prediction for I\ could be obtained by integrating the 
theoretical expression (47) over a very narrow domain near Eq. Clearly, there is no 



way to switch off the radiative corrections in QCD: one has to deal with the pertur- 
bative and nonperturbative contributions simultaneously. The introduction of the 
parameter \x thus becomes mandatory. Eq. (|68| ) then can provide us with one pos- 
sible physical definition of A(/x) (among others) relating this quantity to an integral 
over a physically measurable spectral density. One may rephrase this statement as 
follows. Since quarks are permanently confined the notion of the heavy quark mass 
becomes ambiguous. To eliminate this ambiguity one must explicitly specify the 
procedure of measuring "the heavy quark mass". Any conceivable procedure will 
necessarily involve a cut-off parameter fi much in the same way as the procedure 
defined above, and then A(/z) = Mq — mg(/i). In the "most inclusive" procedure 
when one does not try to separate out any kind of effects, one integrates the tail to 
the kinematical bound E ~ itlq — m q and, therefore, obtains A normalized at the 
scale of energy release. 

Since this question is very important let us look at it from a slightly different 
angle. It had widely been believed that A can be defined as a universal constant. 
The standard definition, being applied to our example, would involve three steps: 
(i) Take the radiative perturbative tail to the left of the shoulder and extrapolate 
it all the way to the point E = E ; (ii) subtract the result from the measured 
spectrum; (iii) integrate the difference over dE with the weig ht function (E-E p Q hys ). 
The elastic peak drops out and the remaining integral is equal to r (t/ 2 /2)A. It is 
quite clear that this procedure cannot be carried out consistently - there exists no 
unambiguous way to extrapolate the perturbative tail too close to Eq Hvs , the end- 
point of the spectrum Q. Our procedure, with the normalization point [i introduced 

7 The standard routine corresponds to using the literal perturbative expression for this tail with 
the non-running strong coupling (for one-loop calculations) , or accounting for the first term in the 
expansion of a s (k) in terms of a s (m) (in the two- loop ones). 
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explicitly, is free from this ambiguity. We will further comment on this issue in Sect. 
5.1 where we discuss the possibility of measuring A(/x) in the inclusive B — > X c lv 
decays. 

In practice, the \i dependence of A(/x) may turn out to be rather weak. This is 
the case if the spectral density is such as shown in Fig. 4, where the contribution 
of the first excitations (lying within ~ Aqcd from Eq VS ) is numerically much larger 
than the radiative tail representing (at least in the sense of duality) high excitations. 
It is quite clear that if the physical spectral density resembles that of Fig. 4 and \i = 
several units x Aqcd, the running A(/x) is rather insensitive to the particular choice 
of \x. As known from QCD sum rules it is just this situation which occurs for the 
standard quark and gluon condensates (the so-called practical version of Wilson's 
OPE). 

Still, even if the /i dependence of A is practically weak, conceptually it is impos- 
sible to define A in the limit \i — > 0. Physically it is quite clear from the discussion 
presented above. This consideration can be thought of as an illustration to a more 
formal argument presented recently |28, f29| . 



If one replaces E in the numerator of Eq. ( |5TD by E (the non-relativistic ap- 
proximation) one arrives at the formula obtained previously in Ref. p4j. Notice 
that the total prediction for 1\ is, of course, \i independent. Eq. ( p7j ) shows how 
A(/i) changes under the variation of the normalization point, 

5A = 5^^. (69) 

9 7T 

The numerical coefficient in front of Sfi is slightly different from that found in Ref. 
(16/9 versus 2-k/3). The reason is obvious - we use here a different procedure 



for defining A(/i) compared to that suggested in PBf : introducing a gluon mass A 
"switches off" the perturbative tail in a 'soft' rather than 'hard' way at fi = A. The 
fact itself of the presence of a linear (in //) renormalization is obviously common for 
all proper procedures. In other words, running of A is not logarithmic, but power- 
like. This fact can be traced back to the mixing between the operators Q(ivD)Q and 
QQ established in [pg| . Numerically S A ~ 0.1 GeV if /i changes from 1 to 1.5 GeV, 
i.e. gK(ji)/K(ji) ~ 0.2. 

By the same token, using the third sum rule, Eq. (§2), one can give a physical 
definition of /i^-(/x), 



r E phys 

= / ° (El hys - Ef^-dE = 
Jo dE 

VM+V / = 2 dE 

JO 671 EqWIq 



(70) 
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where by definition 

/40) 



3r v 



7 phys 



dT 
dE 



^phys _ E y dE 



(71) 



The [i dependence is then obtained as follows: 

4a 



3tt 



V ~ 0.1 GeV 2 if V ~ 1 GeV 



Because the integral of the perturbative tail does not depend on the particular heavy 
flavor hadron in the initial state, this equation is equivalent to the corresponding 
power mixing of the kinetic energy operator with the leading one: 



d 
dfi 2 



Q (vf) 2 Q 



3tt 



QQ 



(72) 



where the numerical value of the coefficient quoted above (4/37r) refers to this par- 
ticular way of introducing the renormalization point. 
The fact that the operators Q (ti) 2 Q 



and 



QQ 



mix with each other, and 



the mixing is quadratic in /i, has been noted some time ago in lattice calculations 
(G. Martinelli raised this question in a private discussion with M.S. in 1991; see 

])• 



3.5 Sum rules for the form factors at zero recoil 

We continue investigating our toy model with the aim of establishing sum rules for 
the form factors at zero recoil. These sum rules will allow us to find corrections to 
the elastic form factor at zero recoil in terms of inelastic contributions. 

The starting idea is to consider the kinematical point q = and go close to 
AM = Mb — Me,. In other words we abandon the case q 2 = and turn to q 2 ^ 
where q is now the momentum transfer carried away by our hypothetical quantum. 
Since q = the ground state B meson cannot decay into a P-wave state. Unlike 
the case of the Bjorken sum rule there is no "external" vector v. The fact that we 
are at zero recoil implies that inelastic contributions are suppressed as A 2 /m 2 , i.e. 
produce effect of the same order of magnitude as the corrections to the elastic form 
factor. 

In the case when q^O the strength of the contribution of the excited resonance 
is proportional to q 2 /M|, ||21|| . In our kinematics q = and the only relevant velocity 
is provided by the primordial motion of the c quark inside D (more exactly, it is the 
'difference' between the quark motions in B and D that counts). The limit q — > 
is obtained, qualitatively, if instead of the "external" velocity we use that of the 
primordial motion, 

qVMl - a^/m 2 . 

8 We hope that this rather clumsy notation will not cause confusion: ^ is the matrix element 
of the kinetic energy operator while jj, (with no subscript) is a normalization point. 
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This analogy nicely illustrates why the residues of the excited resonances are pro- 
portional to Aq CD /M|, and match the picture of Fig. 4 where the contribution from 
the excited resonances is proportional to q 2 /M|,. 

Let us sketch how the sum rules at zero recoil emerge technically. To this end 
we again consider the T product ( pE3D sandwiched between the B meson state. The 
contribution of a given hadronic state in this amplitude is proportional to 

(73) 



(M Q -q ) 2 -M 2 (AM-q )(AM-q + 2M q ) 

It is convenient to introduce the variable 

e = AM - g • (74) 

We are interested in singularities of the amplitude at e M q . The second pole 
in Eq. (|73| ) at e = —2M q is a reflection of a distant singularity corresponding to 
the transition B + D — > <f)* where the asterisk marks the virtual <fi quantum. We 
assume, as usual, the hierarchy A QCD C e < m Q,q- Correspondingly, we will expand 
in AqcdA an d in e/m. In particular, the factor responsible for the second pole in 
eq. ( |T5| ) becomes 

1 . .:. .... . 



AM -q + 2M q 2M g y 2M q 



The second and higher order terms in the brackets can be omitted since they lead to 
a non-singular (at e = 0) expression; the non-singular expressions have no imaginary 
part and are irrelevant for our purposes. 

Next one considers the theoretical expression for the quark level transition op- 
erator, 

(m Q - m q - q ){m Q + m q - q ) + (tt 2 + 2m Q ir - 2q n ) ' ' 

Our task is to expand the transition operator in AqcdA an d in e / m Q,q an d 
then to compare the terms singular in 1/e in this theoretical expansion with the 
phenomenological expression obtained in the language of the resonance saturation. 
A technical point which deserves mentioning right at the beginning is a mismatch 
in the definitions of e. The theoretical expression ([F5|) is phrased in terms of the 
quark mass differences without reference to mass differences of mesons. Since we 
would like to get the sum rules written in terms of the physical excitation energies 
(measured from the mass of the lowest-lying meson state) we have to express Eq. 
([T5|) in terms of AM rather than Am before expanding it. As we will see shortly, 
for our purposes it is necessary to keep all effects through order A 3 /m 3 ; those of 
order A 4 / m 4 and higher can be neglected. 

The expansion for the meson mass in inverse powers of the heavy quark mass 
for spinless quarks considered here is given in Eq. (|33|). It is worth reminding that 
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p 3 denned in Eq. ( p3D is the second order perturbation in 1/rriQ, and as such is 
positive-definite. Using this mass formula one gets 

A A, = A m(l-^ + /^i) (76) 
\ 2m Q m q ^m A Q m l q ) 

where the following short-hand notation is used: 

(tt 2 )o = (H Q \Qn 2 Q\H Q ) mQ=00 . 

Substituting Eq. (|T6|) into Eq. (ff5|) and expanding in A/e and e/rriQ yq using 
Eq. ([35]) we arrive at 

_ {HQl f m = i_L f <g Q> _ <fM) + (77, 

plus terms of higher order in 1/e and in Aqcd- 

This theoretical expression is to be confronted now with what one obtains from 
saturating the amplitude at hand with meson poles. From Eq. (|73|) it is not difficult 
to see that 



- (H Q \f\H Q ) = J*2- £ (iif + + 
where is the form factor for the transition — > HI induced by the vertex 



(7f 



(Hfl ,- Ql ^ = , (79) 

M W is the mass of i-th state ( M<® = M g is the mass of ground state) and 

et = Mf - Mf (80) 

is the excitation energy of the i-th state; all form factors are taken at the zero recoil 
point, where the meson produced in the transition Q — > q is at rest in the rest frame 
ofH Q . 

Comparing Eqs. ([78]) and (|77j ) and using Eq. (^) we find that 

E Ff = l - — !— (H Q 1 2mQQ7r 2 Q\H Q ) (J-+*) + O (±) (81) 
i=ti,... 2M Q \ 2m Q 2m V Vm ' 

and 

y uF 2 = p 3 ( Am ) 2 + ( M . (82 ) 
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Note that the elastic pole gives no contribution in Eq. fl82|). Hence we conclude 
that the residues of the excited states are proportional to A 2 . Transferring then all 
excited states in Eq. (^T|) to the right-hand side we observe that the square of the 
elastic form factor receives corrections of order A 2 , both, of local nature and due to 
excited states. Since all excitation energies and all residues are positive an obvious 
inequality holds, 

^ ^ 9 1 o (Am) 2 . . 

where e\ is the excitation energy of the first excited state. Therefore 

1 (H Q \2rn Q Q^Q\H Q ) (^ + ±-\ - < F 2 < 



2M Hq ' wi ^ ^' V 2m Q 2m V ^e x ra Q m\ 



< 1 - 2^<^ 2m ^ ff2Q i^(w: + «l (84) 



where F is the form factor of the "elastic" transition Hq — > H q . 

It is worth noting that the explicit form of the corrections, in particular to the 
first sum rule ( |SlD and, therefore, the "local" terms oc (n 2 ) in Eq. ([84]) depend on 
the structure of the "weak" current considered and refer to the case of the scalar 
vertex. Should we use the vector current, coefficients in the sum rules would take 
a form leading to Fq — 1 at m q = itlq in accord with the exact conservation of the 
vector current for equal masses. 

Concluding this section let us mention a convenient computational device. It 
is helpful to let the initial quark mass go to infinity and retain corrections only in 
l/m q . In this way one removes nonperturbative corrections originating in the initial 
state. The results referring to finite mq can be simply reconstructed at the very end. 
On the other hand, in the opposite limit, m q ^> mq one suppresses nonperturbative 
effects in the final state; in this way it is convenient to obtain relations for static 
hadronic quantities. 

Similar sum rules at zero recoil in real QCD will be discussed in Sect. 4.2. 



4 Real QCD 

We proceed now to discuss the sum rules emerging in QCD for processes of the type 
B — > XJu. It is clear that the approximation m;, — m c mfc iC is not suitable in 
this case; one still can reach the SV limit, however, by using the fact that q 2 is not 
necessarily zero in this transition (from now on q is the momentum of the lepton 
pair). Indeed, if q 2 is close to its maximal value, 

q 2 max = (M B - M D f , 
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the D meson velocity is small. At the maximal value of q 2 the velocity vanishes. 
The velocity v = —q/E D is related to q 2 as follows 



E D -M, 



D 



M 



D 



(M B - M D ) 2 - q 2 
2M B M D 



J5) 



whereas at the quark level 



1 



(m b - m c ) 2 - q 2 
2mbm c 



(86) 



For a sizable fraction of events measured in the semi-leptonic B decays the values of 
q 2 are such that these events actually do belong to the SV limit (i.e. v is small). An 
indirect proof of the relevance of the SV limit to the inclusive semi-leptonic decays 
of -B's comes from the fact that about 65% of the total semi-leptonic rate is given by 
the "elastic" transitions to D and D* 
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The analysis below is carried out under 
the assumption that the right-hand side in Eq. ( jS5|) is small. Even though we do 
not assume that m;, — m c <C m^.c both quarks, b and c, will be treated as heavy, 
m b , c > Aqcd- 

Needless to say that the proximity of q 2 to g^ aa , can be realized in different ways; 
for instance, one can merely put q = - this is especially convenient if we are 
interested in the zero recoil point - or one can keep q ^ 0, but small and study the 
terms proportional to q 2 . This yields a practically realizable method of measuring 
A(/i) in the semileptonic decays B — > X c ev. We will consider first the simplest sum 
rule for the total decay width analogous to eq. ( |8*ID in Sect. 2.5. Surprisingly, this 
analysis produces a lower bound on the deviation from unity in the B — ► D* elastic 
form factor at zero recoil which does not quite agree with the previous estimates 
obtained by a different method (|0| (see also Ref. ||36|| ). The result is of a paramount 
importance for the experimental determination of |V^&|, the CKM matrix element, 
from the exclusive decay B — > D*ev. Then we turn to an analog of Voloshin's sum 
rule which appears to be a promising tool for extracting A(/x). 

The relation for A as the quantity measuring the mass difference between the 
heavy flavor hadron and the heavy quark has been obtained in Ref. by analyzing 
the heavy quark distribution function appearing in the SV limit for the final state 
quark. It is convenient to rewrite the expressions obtained in Ref. 

Consider, first, Eq. (75) of Ref. [[15 



151 as follows. 



2M 



Im (Hq\T\Hq) 



71 



Hq 



+S'(x) 



Am 2 n A 



5{x)[l- 1 -^ J dyy- 2 G{y)\ + 



'qo(n 2 ) , 1 q 



2m 2 Q A 



+ 



3 m 2 



dy y G(y) 



-2 



G(x) 



where G is the temporal distribution function defined there. The parameter x in 
this expression measures the energy with respect to the quark boundary Eq, rather 
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than to E' Q 



phys 



_ E - E ml-m 2 c + q 2 

X — == , En 



A 2m b 

It is clear then that the 5' term in Eq. (|87|) must merely shift the argument of the 
S term to its physical value, x phys = (E — E$ ys )/A. 

Limiting ourselves to the effects linear in A we get, as a consequence of this 
requirement, 

f ~G(x) = - 3 -. (88) 

J-oo X 2 

This sum rule constraining the temporal distribution function is interesting by itself. 
We can transform it further using its integral representation, Eq. (76) of Ref. fi~5 



G(y) = -i= / dte~ iytA e imbt {B\b(t = 0,x = 0) tt, Te" 1 ^ 1 ''^' tt* &(*,£= 0) \B) 
27rA J 

(89) 

(the factor exp (im^t) accounts for the explicit time dependence associated with the 
rest energy mj, that has not been factored out here). Multiplying this relation by 
y^ 1 and integrating over y we arrive at 

_ = 2 r°° rdt_ -i*e^ |5( } -ifiMW n b u $ = )\ B ) . 

3io ( J 2tt 2M b ^ 1 v ; v ' n 1 

(90) 

Thus A is expressed in terms of a non-local correlator of heavy quark currents (see 
also Ref. [|28|]). Let us remind that the definition (|S6| ) involves local operators made 
of light fields. Its normalization point dependence can easily be traced formally 
through the properties of the path ordered exponent which, being the field operator, 
requires specification of the normalization point for the gauge fields. In Sect. 3.4 we 
have illustrated the renormalization point dependence using the saturation of the 
correlator by intermediate states, which is equivalent to calculating the correlator 
via dispersion relations. 

It will be demonstrated in Sect. 6 that matrix elements of the type entering 
Eq. ( |90|) have simple interpretation in usual quantum mechanics using the first 
quantized language. Here we illustrate it for the case under consideration. Eq. ( |90| ) 
is written in the heavy quark limit where all corrections 1/m are neglected. To this 
accuracy the time ordered exponent is nothing but the (time) correlation function 
of a nonrelativistic heavy quark in the external gluon field: 

(Q(0)Q(x)) A ^ = Te-^oMt')^' 5 3^ e im Q t _ (gi) 
Then the matrix element in Eq. (|90|) can be written as 

j,m h t t 

(B\b{t = 0,x = 0) tt, Te -iJ Mt'm n i b{t,x = Q)\B) = 



c 



2M 



B 



I d 3 x-±-(B\b7r l Q(t = 0,x) Q 7c t b (t, x = 0)\B) ■ e -^" m ^ (92) 
J Z1V1 D 
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where Q is a heavy quark, b or c. This form is convenient for the transition to 
quantum mechanical notations (see Eq. ( |157| ) in Sect. 6) in which this matrix element 
becomes 

(B\iTi(t = 0) 7Tj(i) \B)qm = J2\( B \^\ n )QM\ 2 ^ itnt ■ (93) 

n 

Here 7r is the momentum operator of the heavy quark; matrix elements with the sub- 
script QM are understood in the quantum mechanical sense, namely, only the states 
with zero total momentum are considered and the nonrelativistic normalization is 
assumed. e n are the the excitation energies of the heavy meson, e n ~ M Bn — M B ~ 
Mp n — My). Integrating the correlator over t and £ according to Eq. fl9"0| ) we thus 

HE K T''tr' 2 ■ < 94 » 



The expressions (|90|) , (|94|) are relations that give an alternative to Eq. (pq), 
phenomenological definition of A in terms of measurable correlators. In Eq. (|94| ) the 
sum over n extends up to excitation energy e n = [i and in this way yields A(/z). 

The quantum mechanical derivation of Eq. ( P^ ) as a relation between the mass 
of the heavy quark and the total mass of the quantum mechanical system containing 
heavy quark, will be given in Appendix. 



4.1 Sum rules at zero recoil: generalities 

Our analysis of the B — > X c ev problem at zero recoil will parallel the corresponding 
consideration carried out in the toy model of Sect. 3.5. The presence of spin is 
a technicality which can easily be incorporated. As a matter of fact, all formulae 
necessary for derivation of the first and the second sum rules exist in the literature; 
we will borrow them from Ref. [Hj as well as all relevant notations. 
The point q = represents zero recoil. Then the transition operator 

% v = ij e-^dx T{jI(x)U0)} (95) 

for the b — > c transitions, with j M = cT^b , = 7 M (1 — 75), can be presented (in 
the tree approximation) in the form of the following expansion: 

T KV (u a. 1 ^ 1 V- ^^0 + ^ + ^/2)^ % 



with 



k = m b - g - 

The operator product expansion ( p^) is justified provided that 

Aqcd < \m c -k \. 
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In other words, the expansion (RB) is a series in Aqcd/^c ~ ko). At the same 
time, apart from the poles l/(m c — ko) it obviously contains powers of l/(m c + k ) 
which develop "distant" singularities at k = —m c . We want these singularities 
corresponding to the propagation of the antiquark c to be indeed distant so that the 
dispersion integrals we will be dealing with do not stretch up to these c containing 
states. To this end we must impose a second condition on \m c — k \, namely 

\m c — ko\ <C m c . 

Once this condition is imposed we expand in powers of (m c — ko)/m c and 
Aqcd/(^ c — ^o)- The result is then ordered with respect to the powers of l/(m c — ko). 
The terms non-singular in (m c — ko) are irrelevant and can be discarded. Each 
particular power l/(m c — ko) n+1 in the expansion leads to a sum rule with the 
weight function oc (m c — ko) n . 

Let us sketch the basic elements of the procedure in some detail. We start from 
a series in l/(m c — k ). The next step is averaging of over the B meson state, 

h iaf = -^-{B\t lv \B). (97) 
The hadronic tensor h^ u consists of different kinematical structures j|, 



h^ u = -hg^ v + h 2 v fM v u - ihtf. ia , o #v /ip + h A q^q u + h^v u + q v v^). (98) 

Moreover, the invariant hadronic functions hi to h$ depend on two variables, go and 
q 2 , or q and \q\. For q = only one variable survives, and only two of five tensor 
structures in h^ v are independent. 

Each of these hadronic invariant functions satisfies a dispersion relation in go, 

k(*>) = - f WM ° )d ^ + Polynomial (99) 



2vr J g - % 
where Wi are observable structure functions, 

Wi = 2 Im hi . 

This dispersion representation for hi assumes, as usual, that the integral on 
the right-hand side runs over all cuts that the transition operator may have. The 
general structure of the cuts in the complex g plane is rather sophisticated; the 
issue deserves a special discussion since it is not always properly understood. 

The structure of the cuts of the functions hi(q ) is shown in Fig. 5. The part 
accessible in the decay channel of the B mesons covers the interval [0, Mg — Md]. 
The dispersion integral (p9|) can be written as a sum of two integrals, 



2vr Jo g - go 27r J a g - g 
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where the domain A consists of two subdomains, g < (AJ and q > M B + M D 
(A 2 ). For real decays we are interested only in the first integral since the second one, 
rather than describing the B — > X c ev decay, refers to other physical processes. The 
subdomain A\ actually describes a similar b — > c amplitude, yet with negative go, 
and can be called the lower cut (at go < —Mb — Mr, it contains also the g 2 cut and 
for g < — 3Mb — Md the w-channel contribution is present as well). The integral 
over A 2 , on the other hand, will be referred to as the integral over the distant cuts. 
Two kinds of problems are encountered in evaluating the total dispersion integral. 
The first one emerges due to the fact that for real decay kinematics one has only 
g > 0; therefore, say, for calculating the total width one has not the integral over the 
whole physical cut, but needs to consider the smaller interval without the subdomain 
A\. The corresponding problems of separating the contribution of the same type of 
intermediate states, but at different values of g are usually referred to as "local" 
duality. In the context of the present paper this is however not very important. For 
in our sum rules, from a purely theoretical point of view, it does not matter whether 
a particular transition can be measured in a real experiment or not; e.g. the lower 
bounds we will discuss rely only on the positivity of the corresponding transition 
probabilities. 

There is generally another complication associated with the integral over the 
distant cuts (in particular, subdomain A 2 ) corresponding to quite different interme- 
diate states. The problem of isolating these contributions can be generically referred 
to as "global" duality. 

Both contributions thus represent a contamination for real decays. Fortunately, 
this contamination is irrelevant for our analysis. Indeed, to address the contamina- 
tion due to the "lower" cut, let us choose the "reference" point of go between the 
cuts (see Fig. 5), close to Mb — Mp, 

go = M B - M D -e (101) 

where e is a negative number, 

m D » -e » Aqcd 

When calculating the functions hi at the quark level from the operator product 
expansion we get a similar dispersion relation for the OPE coefficients (with the 
meson masses replaced by the quark ones). We then use duality concepts in identi- 
fying the physical relevance of these cuts. Local duality of QCD means that there 
is a one-to-one correspondence between the part of the OPE coefficients originating 
from the lower cut and the corresponding hadronic contribution in the phenomeno- 
logical (hadronic) representation of hi. The validity of local duality can be verified 
by itself by choosing g in the complex plain close to the particular remote cut, 
the lower one for the case at hand. Therefore, we can systematically discard the 
contribution of that cut simultaneously, in the theoretical expression for hi and in 
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the "phenomenological" saturation. In this way we arrive at the relations 

— r w r ark e n de= — r Wi e n de (102) 
2n Jo 2n Jo 

with Aqcd <C <C Mb, where e is the same variable as defined in Eq. (|101| ), but on 
the cut it is positive. It represents the excitation energy. If, additionally, \i < Mb — 
Md the region of integration in Eq. ( |102j ) lies completely in the physical domain. 
The left hand side of Eq. (|102| ) includes the perturbative corrections as well as 
the power-like nonperturbative terms. The local duality we have invoked to discard 
the contribution from the "lower" cut has an accuracy of the type exp{— Const ■ 
/VAqcd}- 

An analogous analysis can be repeated almost verbatim for the contribution 
of the distant cuts to address the question of "global duality". This duality is 
even more transparent physically and is explicit in all perturbative calculations 
and for calculations in "soft" external fields. On the other hand, in principle, its 
accuracy is generally determined by the similar factor depending on m c , namely 
exp{— Const • m c f Aq C d} , because 2m c determines the distance to the remote cut. It 
is important to remember that it is the ratio m c /AQCD, not mt/AQCD that enters, 
at least at zero recoil. In the real world m c /AQCD numerically is not so large, and 
since the constant in the exponent is unknown one may be afraid of an insufficient 
accuracy of the local duality for D's. At present theory provides us with no clues 
as to the value of the constant in the exponential; the degree of possible violations 
of the duality should be established empirically. With this caveat in mind we still 
believe that heavy quark expansion must work well in B — > X c ev. 



4.2 Sample calculation 

After these more general remarks we return to concrete calculations of the theoretical 
part of the sum rules. To find w1 uark we take the discontinuity of the transition 
operator (|96|), 

1 

-disc 7)^ = 27T&r At (/c 7o + m c + /)x 

f>, (t| . m ,) i*Y^w^- r ,. (103 ) 

^ n! (m c + k ) n+1 

where 8^ n '(x) = (d/dx) n S(x). Using this discontinuity in the left hand side of Eq. 



( |102| ) it is not difficult to calculate all moments of the structure functions in the 
leading approximation. For definiteness we will first consider w A , the first struc- 
ture function (see Eq. fl98f) ) in the transition induced by the axial current. (The 
corresponding "elastic" contribution is B — >• D* .) At zero recoil h^ A is singled out 
merely by considering the spatial components of the axial current, 

hf A = hf k A - (1 = 1,2,3). 
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All other structure functions and other currents can be treated in a similar manner. 

Technically it is convenient to carry out the computation in two steps: first one 
obtains auxiliary moments convoluted with the quark value of e, 



m b -m c - q 



c — m c , 



(104) 



and then, at the second stage, these results are converted into the true e moments. 
Notice that in the case at hand e should be defined as 



e = M B - M D * - q 

since the lowest-lying state produced is D*, not D. Then, according to Eqs. fl26[)-(|29"D 



(M B - m b ) - (M D * - m c ) 



I 2 2 \ ( 1 



\2m r 2mt 



3m c 



/4 + 



'A 4 



(105) 

where /j^ and /i^ defined in Eq. (|l|) are the asymptotic expectation values of the 
kinetic and chromomagnetic operators. In the estimates of the neglected terms m 
generically denotes both m c and m&. The mass shift 5a determines the difference 
in the threshold energy between the real hadrons at zero recoil and the quark mass 
difference, i.e. it is a direct zero recoil analog of A in Eq. fl53|). 
Let us define moments of the structure functions Wj as 



J« = - / e n Wi (e) de . (106) 

Z7T J 

In the leading non-trivial approximation we get for the moments of W\ 
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a k {f + 3B)( 



1 



2m r 2rrih 



2m, 



[aB,a k ] 



b) + 0(AS/ 



m 



(107) 



where B denotes the chromomagnetic field, and the last equation refers to n > 2. It is 
worth noting that to this order in 1/rriQ only chromomagnetic field appears explicitly 
for n > 2 although in the original expansion ( 103 ) the field tensor G contains both 
chromomagnetic and chromoelectric components; the chromoelectric field cancels 
out. The underlying reason for that will become apparent shortly (see Sect. 6). 
Averaging (...) over the initial hadron state is understood as (2Mb)~ 1 (B\...\B) in 
all expressions in the right hand side. In the derivation of the above predictions for 
the moments we used the QCD equations of motion (TT3]). 

The structure of the solution of Eqs. (|107|) for the excitation function W\{e) 
is quite transparent. The first equation tells us that the sum of all probabilities 
is equal to unity up to small corrections ~ (9(AQ CD /m 2 ). On the other hand, all 
higher moments I n explicitly start with terms of order Aq CD • (AqcdA^) 2 - Since the 
scale for the excitation energies e« is given [] by Aqcd one immediately concludes 
that the probabilities of transitions to the excited states all scale like Aq CD /m 2 . To 
saturate the first sum rule one then needs state (s) which do not contribute to the 
higher moments and are produced with practically unit probability; the only way 
to satisfy this constraint is to saturate by the final states D and D* with the masses 



up to corrections vanishing in the limit mb,m c — > oo. Moreover, these "elastic" 
transition amplitudes must be equal to unity up to terms inversely proportional to 
the square of the heavy quark masses - a fact observed originally in Ref. |23| and 
now known as Luke's theorem |42|. Although the expressions above are derived for 
the axial current, similar results are valid for the vector currents as well, where the 
lowest-lying final state is D, not D*. In this way one obtains also the statement of 
the heavy flavor symmetry in the spectrum of hadrons. 

The second of Eqs. ( |107| ) corresponding to n = 1 has a special status. When 
written in terms of the excitation energies e, counted from the quark threshold, 
the right-hand side does contain terms of order 0(A^ CD /m) which are expressed in 
terms of /x 2 and /i^; on the other hand the structure of the solution described above 
implies that the contribution of the excited states in the phenomenological part is 
only of the order of AQ CD /m 2 . Therefore, this leading term is to be completely 
saturated by the elastic peak that resides not at e = but is rather shifted by 
the amount 5a = e — £■ This condition unambiguously determines the leading, 
0(AQ CD /m), correction to the masses of the heavy flavor hadrons which has been 



anticipated in Eq. 



Similarly, to order 0(h% CD /m) the sum rule for the second 



Effectively the same refers even to the thresholds associated with the £)W + pion(s) in the 
chiral limit which, strictly speaking, have no excitation gap. This is true due to the fact that the 
corresponding amplitudes are proportional to the pion momentum; they can produce only chiral 
logarithms and do not change powers of mass in the analysis, see Ref. fH[. 
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moment yields the "local" 1/m 2 term in the effective Hamiltonian, Eq. and 
expresses the non-local correlators p 3 as the sum over inelastic probabilities. 

It is not difficult to obtain a general expression for the function w AA as it emerges 
from its moments I^ AA . The inelastic part of w(e) appears at the 1/m 2 level and 
is given by the Fourier transform of the time-dependent correlation functions of the 
leading operators in the effective Lagrangian: 



e w 



AAi 

1 I 



3 



X e -» IoMt',o)dt' 



2m c 

1 lxl 



a k (n 2 + aB)(-^- - -^-) - -^—[a k , SB] 
2m c zmf, 2m c 



b(t,x = 0)\B)' . 

(108) 

The prime here means subtraction of the "elastic" contribution (see the expression 
for I 2 and the discussion below). 

Equation ( |108| ) has a very transparent quantum-mechanical meaning. It corre- 
sponds to the following picture. At time t — the Hamiltonian of the system under 
consideration is suddenly changed by adding a perturbation 



/ 1 



\2m c 2mb 



(ir 2 + aB) . (109) 



Simultaneously the original wave function is changed due to the spin flip. This effect 
is represented by the term (— l/2m c ) [o^, aB] where k marks the spatial component of 
the axial current. Had we considered the time component of the vector currents the 
latter term would be absent, and the entire perturbation would reduce to Eq. ( |l(Jij| ). 
At time t the perturbation is switched off. The second order term in this combined 
perturbation yields us the excitation probabilities. This interpretation brings us 
very close to Lipkin's quantum-mechanical formalism which will be discussed in 
some detail in Sect. 6. 

Transitions induced by the time component of the vector current are given by 
the following combination of the structure functions: 

w w = _ w vv + w vv + q 2 w vv + %w vv (no) 
The moments of w vv defined in analogy to Eq. (|106|) look even simpler: 

tr = a {(*■ + sB) (± - ^) + (-(Be) +2a • e x ?) ) ( ' ) } b) + 



+ Sv + 0(A l QCD / 



m 



36 



- S 2 V + 0(A^ CD /m 3 



rVV 



,7T 



BB){ 



2m c 2111b ' 



X 



,7T 



b) 



0(A n Q +? D /m 3 ) 



'III) 



"2m c 2171b' 

Note that for the vector current the lowest charmed state is the D meson and, 
therefore, e is to be defined now as e = Mb — Mb — q ; the value of 5y then differs 
from 5a, arid the second sum rule (the expression for the first moment) yields now 



(M B - ,n t ) - (M D - ™J = -« - &) (-L - -L) 



+ 



p 3 



pIg 



Pi 



P\ 



4m? 



Ami 



+ 0(A^ CD /m 3 ) (112) 



instead of Eg. (|T05|) . 

It is not difficult to trace a straightforward parallel between the derivation given 
above, Eq. ( |108|) , and the analysis of the effects of the heavy quark motion in the SV 
limit carried out in Ref. |[15|] that led us to the temporal distribution function. In that 
case, however, inelastic excitations appeared proportional to the velocity squared of 
the final quark, and it was the correlator of the spatial momentum operators tt that 
emerged. 

A remark is in order here. Obviously, one can, in principle, calculate higher 
order power corrections to the sum rules ( |107| ), ( |1 1 1|) . It is possible to see that - 
in what concerns the inelastic contributions to the structure function - w(e) will 
order by order reproduce the successive terms of the ordinary quantum-mechanical 
perturbation theory corresponding to both next order iterations of the leading terms 
as well as to perturbations representing subleading power operators in the effective 
heavy quark Lagrangian; the local relativistic corrections in currents appear as well. 
Similarly, one can consider the sum rules in the case of nonzero recoil. This cor- 
respondence to the perturbation theory in quantum mechanics will be discussed in 
Sect. 6. 

To conclude this section, let us say a few words elucidating why this general 
analysis based on the sum rules is important although it might seem to lead to 
no new results beyond the picture of the standard perturbation theory in quantum 
mechanics. 

First, we shall see in the subsequent sections that this approach allows one to 
obtain useful bounds on the transition amplitudes and on the 'static' matrix elements 
governing nonperturbative corrections. 
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Second, it demonstrates in a transparent and unambiguous way the necessity of 
introducing the infrared normalization point in addressing nonperturbative effects 
and clarifies both the physical meaning and the qualitative trend of the \x depen- 
dence. This is important in view of the existing misinterpretation of this problem 
in applications of HQET. 

The last, but not the least aspect: we clearly see here that such assumptions 
about QCD as the validity of global duality discussed above, that sometimes are 
naively believed to be specific only for the OPE-based approach to inclusive transi- 
tions, are actually the most general requirement inherent to any consistent consider- 
ation. If for some particular physical reason the quark-hadron correspondence used 
above were modified - for example there were a sizable "leak" in dispersion integrals 
from distant cuts, it would immediately lead to new contributions in the right-hand 
side of the sum rules; this would result in physically observable corrections not 
contained in the expansions that can be obtained in HQET. 

4.3 A bound on the form factor at zero recoil from the sum 
rules 

The general theory developed above is applied in this section to derive a lower bound 
on the deviation of the "elastic" form factor at zero recoil from unity. To this end 
we analyze the first sum rule. 

We consider, for definiteness, transitions generated by the axial current, 

= 07^756 , 

i.e. the transitions of the type B — > D* and B — > excitations of the vector mesons. 
Practically they are most important in the exclusive approach to a determination 
of \Vcb\. These transitions are induced by the axial- vector current and as was 
mentioned above it is most convenient to focus on the spatial component of this 
current. For the spatial components of the current only hi survives. 

To get the first sum rule at zero recoil we use the first of Eqs. ( |107[) which to 
order 1/m 2 reads as 

1/*^) = 1-14-^(^ + ^+2) . (113) 

The chromomagnetic parameter is known experimentally: 

/4^(M^-il4) = 0.37GeV 2 . 
The sum rule stemming from Eq. (|113| ) obviously takes the form 

r _B->D* T Z_^i ^B^excitations 



3 m 2 4 ym 2 nit 3m c mb J 

where the sum on the left-hand side runs over all excited states and all form factors 
are taken at the zero recoil point. This is a perfect analog of eq. (|8ip . The form 
factor B — > D* at zero recoil is defined as 



(B\A k \D m ) = y AMbMd* 

where is the polarization of D* meson. 

If all terms (9(Aq CD ) are switched off higher states cannot be excited at zero 
recoil - only the elastic B — > D* transition survives - and we arrive at the well- 
known result that 

Fb-*d* — 1 (zero recoil) , 

the statement of the heavy quark (or Isgur-Wise symmetry first noted in the 
SV limit in Ref. |23j (see also |44]| ). Including C(Aq CD ) terms we start exciting 



higher states; all transition form factors squared are proportional to AQ CD /m 2 . Si- 
multaneously the form factor of the elastic transition shifts from unity. 

Both power corrections on the right-hand side are negative. What is crucial is 
the fact that the contribution of the excited states is strictly positive. Transferring 
them to the right-hand side we arrive at the following lower bound 



F 2 ^ , > + ti Ok ( _L _|_ _L + f ) . (H5) 

L> ' L> 3m 2 4 \m 2 ml 3m c mf }/ 



Since /i 2 > p? G (see below) we find that the (absolute value of) the deviation of 
the elastic form factor Fb-,d* from unity at zero recoil is definitely larger than 

ML - Ml 

B Q 9 B ~ 0.035 . 

As was mentioned, the second term on the right-hand side of eq. ( |115| ) is also 
positive, /i 2 > /Mq. We will rederive this inequality within the framework of sum 



rules themselves in the next section. Previously it was obtained in this form in [45 



where the quantum-mechanical argument of Ref. |15j was extended. The quantum- 
mechanical line of reasoning is applicable at a low normalization point. 

The inequality n\ > p? G is in perfect agreement with the most refined QCD sum 



rule calculation of ^ [46[ which lead to /i 2 = 0.6 ± 0.1 GeV 2 ; the recent updated 



analysis along the same line of calculations yielded (see Ref. |T7j ) 

/i 2 = 0.5 ±0.1 GeV 2 . 



If this estimate is accepted then the second term in eq. (|115|) amounts to ~ 1/2 of 
the first one, and the lower bound for the deviation becomes 1 — F^ D , > 0.1. 
The actual deviation is probably twice as large. First, the sum rule derived above 
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neglects perturbative a s corrections. The first order correction to the elastic form 
factor was calculated in Ref. [23|. If in zeroth order in a s the 67^750 axial- vector 



vertex at zero recoil is unity, the first order correction renormalizes it to 

a s fm b + m c m b 8 

Va = 1 H In 

7r \rrii, — m c m c 6 

Numerically t]a — 0.97 if one uses a s normalized at the point p = ^Jm c m h (for 
the recent discussion see Ref. |[48|| ). For the axial- vector current the perturbative 



correction is negative, so that unity in eq. (|1 14|) is replaced by approximately ™ 

'114) is strictly positive, and 



0.95. Then, the contribution of the excited states in eq. 
this also reduces This contribution may be as large as, roughly, the power 

correction on the right-hand side. An estimate of the excited state contribution 

ill where a more detailed numerical 



supporting this statement is given in Ref. 
discussion of all corrections is given. Notice that in our approach the excited state 
contribution replaces a non-local contribution of Ref. . 



We conclude that 1 



>D* 



is definitely larger than 0.1, somewhat beyond the 



window obtained in Ref. [ 4"D|] . The phenomenological impact of this observation is 
discussed in Ref. [|4[ . 

Let us note that in order to use the second sum rule similar to Eq. (|84"D we would 



need to know 0(Aq CD ) terms both in the transition operator and in the relation 
between AM and Am. The corresponding expressions are provided by Eq. (|105|) 



and the second of Eqs. ( 107 ); however the relevant hadronic parameters p^, p\ s and 
non-local correlators p^^G s a are no ^ known yet. 

In a very similar way one can obtain the bound and estimate for the vector form 
factor of the B — > D transition Fb->d at zero recoil. Here only the timelike com- 
ponent of the current contributes, and for this reason the full semileptonic decay 
amplitude is proportional to the lepton masses. Therefore this mode is not advan- 
tageous; the corresponding form factor is measurable (in principle) at zero recoil in 
the B — > D + tv t decays. Taking Tn = T u = 7 we obtain for this case the sum rule 



+ 



E 

=1,2,. 



Fl 



texcit 



1 - 



m c 



--) 



(116) 



Perturbative corrections also differ and now look as follows [23] : 

1 ( m h + m c ln ^fe _ 2 

71 \m b — m c m c 



117) 



The corrections in the case of the vector form factor obviously vanish at m b = 
m c , as they have to in view of the exact conservation of the current in this limit. 
Numerically therefore they are expected to be smaller for vector transitions than for 
axial ones. 



10 For numerical estimate w e use somewhat larger number than the literal value of i]\ above, see 
discussion in Sect. 7 and Eq. (|l89| ). 
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It is worth mentioning that the excitation probabilities entering the sum rules 
( |114j) and ( 116|) are generated separately by the axial or the vector current, respec- 



tively, but not by the V — A current that directly produces the experimental widths. 
Actually at zero velocity transfer the axial and vector currents cannot interfere. 
Therefore for the V — A semileptonic transitions into massless leptons one has just 
to add to Eq. (|114| ) (assuming that no final state identification is attempted) the 
contribution of the cr/jfr current. The sum rule for this current is obtained in the 
next subsection. Combining the two sum rules one gets 

H^ = l -£zJi>-[ de d^l (118) 

ivn^VTLc Je>M D *-M D 2,11 

where the last term representing the inelastic contribution is expressed via the dif- 
ferential semileptonic width at zero recoil: 



wY~ A (e) 8tt 3 1 d 2 Y 



2tt G 2 F \V cb \ 2 q 2 \q\dq 2 dq 



s7 - 1 (119) 



9=0, 10 = M B- M D- 



(q is the momentum of the lepton pair in the process). Equation ( |118| ) is much less 
useful as an upper bound because the main part of the correction to the elastic form 
factor must now come from the presently unknown contribution of excitations. 



4.4 Lower bound on fi\ 

Reversing the line of reasoning used to derive the upper bound on the form factors 
(the lower bound on 1 — F) , we can exploit the very same idea to get a constraint 
on the matrix elements li g and /a 2 . At zero recoil there are only two independent 
structure functions for the correlator of the V — A currents; similar functions can 
be introduced for other weak vertices as well. Choosing a particular current one 
projects out a certain combinations of the structure functions. It is important that 
for the Hermitean conjugated currents in the two weak vertices one gets a definitely 
positive structure function expressed as a sum over certain transition probabilities. 
If an appropriate channel is found where the elastic peak is kinematically absent, the 
theoretical side of the corresponding sum rule will contain no unity term, and start 
with the leading 1/m 2 corrections given just by a linear combination of [i 2 G and li 2 . 
The "phenomenological" side, a sum over the excited states, is positive-definite. In 
this way we arrive at a constraint on the linear combination of li g and fi 2 at hand. 

It is not difficult to find a specific example. Indeed, let us consider the vec- 
tor current, q^^Q, and in particular its space-like components. It is obvious that 
— (l/3)/i^ v = — h\ v is populated only by the excited states - if the initial state Hq 
is the ground state pseudoscalar, the final one, H* must be the axial- vector meson, 
non-degenerate with Hq in the symmetry limit. (In the quark model language one 
would say that H* is a P-wave state.) Using the results of Ref.|lO| (Eq. (Al) ) we 
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find 



h\ V {qo) = - 
e 




4 V m Q m q 3mQm q J 3 



(120) 



"q J 

The expression in the square brackets is equal to the sum over excitations: 

,2 ,,2 



tit-HQ ( 1 I 1 2 ) I 1 = y-^2 fl21 x 

4 3m Q mJ 3 m 2 q Y Ho ~"" 
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and, hence, is always positive, for any values of mg and m q . Being interested in 
the "static" properties of the initial state only, it is convenient, according to the 
comment in the end of Sect. 3.5, to consider the theoretical limit m q ^> vtlq where 
only initial state effects survive []]. Requiring the positivity of Eq. ( |120| ) at m q ^> mg 
we conclude that 

fil-fi 2 G = 4m Q {£|F„ Q ^,| 2 ) (122) 



/4 > /4 ■ (123) 



and, therefore, 



This is literally the same inequality that has been obtained previously [H5J in the 
quantum-mechanical language along the lines suggested in Ref . [jn| . The argument 
presented above can be viewed as a consistent and transparent field-theoretic rein- 
carnation. To be fully consistent, according to the general discussion of Sect. 3.4, 
one must introduce the cutoff in the "phenomenological" integral over the decay 
probabilities from the upper side of excitation energies e. It is most important that 
the integral in the right-hand side of Eq. (|122 ) 



is positive for any normalization point /i, and therefore ensures the validity of the 
inequality (|123|) for operators normalized at arbitrary values of /i, provided that 



the normalization point is consistently introduced in this particular way. For high 
enough \i the contribution of the excited states is given by perturbative expressions 
and the \i dependence is explicitly calculable; of course, for very large \i the inequality 
becomes trivial. We will discuss the issue in more detail in Sect. 7. 

It is instructive to note that the zero recoil sum rule for the x Vk transitions 
(or similar ones where there is no elastic peak) provides us with the direct way of 
determining the evolution of the kinetic energy operator. One proceeds here along 
the same line of reasoning as has been outlined in Sect. 3.4 where we considered the 
small velocity kinematics with q^O. In this case the elastic peak identically vanishes 
for purely kinematical reasons. Therefore the only possible impact of introducing 



n The very same bound can be obtained by considering, say, the correlator of two 175 currents. 
In this case one gets directly the difference fi^ — fj,Q with the coefficient (l/m c + 1/mf,) 2 . 
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the infrared normalization point (say, via the gluon "mass" ) can emerge from the 
gluon emission probability. 

Technically, we can consider the relation ( |121| ) in the perturbation theory. Then 
the hadron state Hq must be replaced by the heavy quark Q, and the excited states 
are q + gluon states. For such initial state the expectation value of operator oG 
vanishes and the sum rule (|121 ) converts into the relation for the perturbative part 



of/4: 

(/4) pert / 1 1 2 \ 1 r d 3 k 1 2 

4 \w?QW? q 3m q m Q ) ' 3^ L<„ 2^(2vr)3 4m q m Q mmkW)l ' 

(124) 

the sum over gluon polarizations is implied. Following the procedure used through- 
out this paper, we have introduced the renormalization point // as the upper limit 
for excitation energy. The calculation of the amplitude in the right hand side is very 
simple: 



(qg\rfQ\Q) = g s qt a 



1 1 \ 1 



e a I + I - i[e" x a 



2m q Ivtlq J \2m g 2mc 



Q . 



(125) 

Taking the square of this amplitude and summing over the gluon polarizations e a we 
get the same dependence on the masses m q , ttiq as in the left hand side of Eq. ( |121| ), 
which is expected on general grounds. In this way we obtain 



(/4 



pert 37T 



U 'V . (126) 



This result can be rewritten as 



d ; Q( l DfQ = ^^QQ (127) 



djJ? 37T 



which coincides with Eq. ([72]) . It is easy to check that the same evolution law 
is obtained for any suitable current and even for a case when heavy quarks were 
spinless. 

Throughout this paper we have phrased our discussion of real QCD in terms of 
transitions where initial states were heavy flavor mesons, namely B. It is clear that 
exactly the same reasoning can be applied for the transitions of heavy baryons, for 
example when the initial hadron is At. The matrix elements are different, of course. 
In particular, the expectation value of the chromomagnetic operator vanishes for 
Moreover, contrary to the meson case no non-trivial lower bound on the kinetic term 
emerges. Therefore it is natural to expect smaller deviations from the symmetry 
limit for both vector and axial-vector form factors in baryons than for mesons. 

An analysis of corrections to Fb-,d* resembling ours in spirit, but not technically, 
has been carried out recently in Ref. [|3(J. There 1 — Fb^d* is expressed in terms 
of some local and non-local expectations values; the latter are unknown. In our 
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analysis the role of the non-local expectation values is played by the contribution 
of the excited states. What is crucial, this contribution is always positive. In |36| 



1 — Fb->d* is found to be positive (good!), and a numerical estimate is presented 
relating 1 — Fb^d* to [i 2 , a parameter that is somewhat more uncertain than \j? G . 
As a matter of fact, the numerical values of /j, 2 accepted in |3(| under the influence 
of some recent claims are in contradiction with the inequality (|123|). 



5 Sum rules at q ^ 



We leave now the point of zero recoil and discuss the sum rules in general situation. 
The first obvious complication is that there are more independent structure functions 
that enter separately the decay rate for any particular current, and each depends on 
two variables for which we choose qo and q 2 . Inelastic processes corresponding to the 
transition to the states other than D and D* are now not necessarily suppressed by 
powers of K/tuq. The n-th moments of the structure functions are therefore propor- 
tional to A . If nonperturbative effects are calculated explicitly through 0(Aq CB ) 
terms, we obtain non-trivial corrections for the first three moments with n — 0, 1 



and 2. We use below some results of the forthcoming paper |n| where more details 
will be given. 

At q 7^ the variable e is defined as 12 



M. 



B 



+ Q' 



qo 



(128) 



and moments of the structure functions are 



1 

2^ 



J dee n Wi(e, q- 



(129) 



where i labels the structure function; they can be considered separately for axial 
current (AA), vector current (VV) and for the interference of the two (AV). 
(A1-A6) obtained in Ref. 



Using Eqs. 



T0|| and expanding the corresponding 
hadronic invariant functions hi in powers of 1/e through terms 1/e 3 one arrives at 
the set of relations which are valid up to terms 0(Aq CB ). We present here these 
relations only for the structure functions u>i,2,3 which contribute to semileptonic 
decays with massless leptons (the most general case will be considered in Ref. |5T 
We have the following sum rules for the zeroth moments: 
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(130) 



12 When vector current is considered, the lowest state appearing is D; we still use the D* energy 
as a reference point keeping in mind that explicit account for the D contribution in the phenomeno- 
logical part of the first moments l\ v is necessary. The fact that the D contribution to semileptonic 
width vanishes at q = for massless leptons justifies such a choice. 
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T (2)AA 
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T (3)AV 
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(131) 
(132) 



Expressions for the VV functions are obtained from the axial ones by replacing 
m c — > —m c ; the structure functions w^' 2 ^ AV and w^ AA ' vv vanish. 

The above equations are analogs of the Bjorken sum rule; however they incor- 
porate nonperturbative effects which appear at l/m^ level; the corrections are not 
universal and differ explicitly for different currents and structure functions. 

The first moments look as follows: 
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(135) 



At q = these relations determine 1/m terms in the masses of heavy mesons. 
Their derivatives with respect to q 2 near zero recoil give the Voloshin's "optical" 
sum rule. Here they are obtained with better accuracy for arbitrary, not necessary 
small, velocity and incorporate I/tuq relative corrections. The latter appear to be 
quite sizable when q is not particularly large. 

The third sum rules, which are relations for the second moments of the structure 
functions, are calculated only in the leading non-trivial approximation. They look 
rather simple and manifestly satisfy the heavy quark symmetry relation |H] : 
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(136) 



All higher moments vanish in our approximation. We used above the notation Ed* 
for the energy of D* and E c for the energy of the c quark in the free quark decay: 



E 



D* 



M 2 D , 



E P 



m 2 + q- 



(137) 
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The quantity 



(M B - m b ) - (E D * - E c ) 



EJ ^ w \2E C 2mJ 3E C 2E C \ E 2 C ) U 2 , 

(138) 

that enters first and second moments, is similar to 5a we have discussed in the case 
of zero recoil; at q ^ however, it is of the order of Aqcd- 

It is important for the analysis of the inclusive widths that the nonperturbative 
expansion of the invariant hadronic functions hi, and, therefore, of the moments of 
the structure functions always runs in the inverse powers of E c rather than m c ; it 
is correlated with the fact that non-physical singularities in go representing distant 
cuts are also separated from the physical cut by ~ E c . It implies that one can use the 
same expansion literally even if m c — > as long as \q\ 3> Aqcd- In turn, this means 
that at E c ~ m& dominating the inclusive width for ^> m c , nonperturbative 
corrections scale like the inverse square of the mass of the heavier decaying quark 

BEB- 

Using the same fact one may hope to decrease l/m c corrections to the determi- 
nation of A and /i 2 from the second and third sum rules, which are rather significant 
numerically in the SV limit, considering decay processes with large recoil. A more 
detailed analysis of the sum rules at non-zero recoil is given in Ref. [Q . 

5.1 The second sum rule at q ^ 0; measuring A(/i) 

We discuss now a few useful applications of the sum rules presented above in the 
the simpler case of the SV kinematics when we can expand moments in g 2 /m 2 . 

If at q — the second sum rule requires the knowledge of (9(Aq CD ) terms, at 
q 7^ (i.e. Aq CD <C \q \ <C M D ) a non-trivial prediction, an analog of Voloshin's sum 
rule p4] , arises in order Aqcd- Higher order corrections are written in Eqs. ( |133|) - 



( |135| ); they will be discussed elsewhere |p(| . We need to consider the average value 
of 5omax — qo- If the value of g 2 , rather than <f 2 , is fixed, this quantity is related 
to Mx c , the average invariant mass squared of the hadronic system produced 0- 
Indeed, then 

M 2 _ m 2 _i_ Q 2 

M 2 Xc = M 2 D + 2M B (q 0max -q ) , g max = B on / - (139) 

If q 2 is fixed then 

M 2 c = M 2 D + 2E c (g 0max - g ) + (g 0ma x - go) 2 , <?omax = M B - J M% + g 2 ; (140) 



13 The corresponding analysis of the average invariant hadronic mass presented in Ref. Q was 
incorrect, see Ref. p| . The average invariant mass of the final state hadrons is not given directly 
by matrix elements of local heavy quark operators at the level of nonperturbative corrections, 
contrary to claims in Ref. B . 



46 



in this case the second moment also contributes, but only at the level of 0(Aq CD ). 
We will use the energy of D* in what follows, and therefore need to define the 
corresponding threshold energy gg max : 



= M B - y/M*, + q* . (141) 

The basic idea is the same as that demonstrated in Sect. 3.4 in the toy model: 
to order Aqcd i n the SV limit the weighted integral over the excited states is pro- 
portional to A(/x). In the previous section we considered the point of zero recoil; 
now we have to shift from this point and consider terms proportional to the square 
of the c quark velocity. The SV limit will be ensured by choosing \q\ <C Mb- 

Let us assume first that all structure functions in eq. (|98| ) are known separately. 
Then it is most convenient to consider the function h\ for axial current transitions. 
If we are aiming at effects linear in A QCD , all fi 2 and [Mq corrections can be neglected; 
we can include them, however. Using Eq. (|133| ) and the similar one for the vector 



current as well as the expansion (|138[) we get 



+ O (g, A^co) (142) 

where the hadronic tensor considered is that induced by the V — A current. In this 
equation q is supposed to be fixed and small as compared to m 2 . 

Let us note the explicit presence of the normalization point \i in the lower limit 
of integration. Contrary to the naive quantum mechanical description, in real QCD 
the structure functions do not vanish when excitation energy becomes larger than a 
typical hadronic mass scale, but rather contain a long tail associated with the gluon 
emission. To reiterate the conclusion of Sect. 3.5: we introduce a normalization 
point /i in such a way that all frequencies smaller than /x can be considered as "soft" 
or inherent to the bound state wave function; at the same time a s (/j,) has to be 
sufficiently small for the perturbative expansion in a s (/x)/7r to make sense. We then 
draw a line at q = Qomax — f 1 (the picture is similar to that of Fig. 4, with dT/dE^ 
replaced by wX~ and by go)- The integral ( |142| ) taken over the range go max — I 1 
to Qomax represents (v 2 /2)A(fi) modulo corrections of higher order in v and in Aqcd- 
The running mass is then defined as rrib{n) = Mb — A(/i). 

Practically it may be not so easy to separate different structure functions from 
each other, which would require the triple decay distribution over q , q 2 and Eg_. A 
similar prediction can be given for double differential distribution in the semileptonic 
decay when the integral over the energy of lepton is considered. The explicit form 
of the sum rule depends again on whether one fixes q 2 or q 2 in the process. Below 
we assume that q 2 is kept fixed. 
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Using the relation 



d 2 T 
dq dq'- 



cb\ 



|2 F 

16vr 4 



/ 2 y 



(143) 



and Eqs. (p3|) , (p4|) , ( ggg ) and ([T38|) one arrives at 



OWOmax 



Qo) 



d 2 T 
dq dq- 



G\ 



\V cb \ 2 {m b - m c ) 



2 m 



2M 2 D 




in/, — in,. 3(m b — ///, )-' J ' Q CD ) ^ ^ 

Again, one can consider radiative corrections. For simplicity we will neglect all 
terms that are ~ Aq CD and higher (then the difference between ggmax an d 9omax 
becomes unimportant), and limit ourselves only by terms ~ q 2 . Including the 
radiative tail we have 



8tt 



G 2 F \V cb \ 2 



(/On 



q\ 



dq (q { 



Omax 



Qo) 



1 d 2 T 
q \ dq^dq 2 



57T 



G 2 F \V cb \ 2 



d 2 Y 



'<?Omax— A* 



\q | dq^dq- 



+ 



<?Omax-M 1 d 2 Y 

dqo{qo max — qo) t^t -. — t^t 
q\ \q\ dqodq 1 



(145) 

The first term in the right hand side, up to the known factor, has the meaning of 
the running value of A if /i is much smaller than \q\: 



8tt 3 



G 2 \V cb \ 2 



'/On 



dg (<?0max - qo) 



d 2 T 



'/On 



\q\ dq dq [ 



[m b — m c 



,2 19 I 



2m 2 



A(ju) (146) 



If 9omax — 9o is still sufficiently large the integrand coincides with the result 
obtained in the perturbative calculation. To the first order in a s it is given by 
the probability to emit gluon with energy uj w gomax — qo- It looks particularly 
simple when gomax — qo ^ \q\- As compared to the case of zero recoil we have 
discussed previously, now the dipole gluon emission appears which is proportional 
to v 2 , and for small uj one can neglect 1/m suppressed amplitude considered before. 



In the SV limit, when q 2 -C m 2 



the results for the radiative corrections obtained 



in Section 3.4 for the toy model are directly applicable to the semileptonic decays. 
Indeed, the gluon can be emitted (absorbed) either by the color charge of the c 
quark or by its magnetic moment. Moreover, there is no interference - the gluon 
emitted by the magnetic moment has to be absorbed by the magnetic moment. As 
long as we consider gluons with momenta much less than q, we can disregard the 
gluon interaction with the magnetic moment. Then we are left with the charge 
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interaction only which is the same for spin-0 bosons and spin-1/2 fermions. As 
a result, the expression for the radiative correction obtained previously in the toy 
model is modified in a minimal way, only due to a slightly different kinematics, and 
for gomax - <?o < \q I we get 

8r ' 1 £r >" - (m„ - m c f - , 16a ' ; . (147) 



G 2 F \V cb \ 2 \q\dq dq 2 ' c 2m 2 c 97r(g 0max - q ) 

The above equation shows how the the value of A obtained from the sum rules 
depends on the parameter /i; it coincides with Eq. (|69|). If one does not introduce 
this explicit cutoff, the value of A would generally scale like a s ■ tuq. 

Thus, the sum rules ( |142| ) and ( |146| ), ( |147| ) can be used to elucidate what is 
actually meant by the heavy quark mass. This question is rather subtle since the 
heavy quark mass is a purely theoretical parameter which is not directly measurable. 
On the other hand, it is a very important parameter, crucial in a wide range of 
questions. 

The sum rules ( |142|) or ( |146| ) express A in terms of the integral over the physically 



observable quantities. Therefore, we have a suitable phenomenological definition of 
A and, through this quantity, the heavy quark mass. Of course, both of them depend 
explicitly on the normalization point. 

Although the above equations yield A(/i) in terms of the excitation distributions 
which are, in principle, experimentally measurable, this does not mean that it can 
be easily measured in practice. Needless to say that it has not been measured so 



far. Still, the expression (|142|) , combined with the Bjorken sum rule, implies [51] a 



lower bound on A(/i) (see also P4]| ) which turns out to be quite restrictive, 

A>2A 1 (p 2 -i), (148) 

where Ai is the mass difference between the first excitation of D and the D meson; 
p 2 is the slope of the Isgur-Wise function [^0, 21 1. Numerically the right-hand side 



of Eq. ( |148| ) is close to 0.5 GeV. Further details can be found in Ref. [5T| . 

5.2 The third sum rule in the SV kinematics 

In a similar manner one can use the third sum rule to relate the kinetic operator 
to the average value of (qo ma x ~ go) 2 - The value of p 2 , can be extracted in a model- 
independent way from the sum rule similar to Eq. ([H]) if the double differential 
measurements are used; say, for small velocity events 

Hi (ji) = 3r~V 2 / — — (q 0max - g ) 2 dq (149) 

J io ma x-i j - «go«g 



14 Virtual corrections lead to the term ~ <5(qomax — <7o) which vanishes in the moment we consider, 
but for the zcrot h moment in the first sum rule it cancels the logarithmic singularity following from 



expression (147) 
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where v = \q\/m c or \q\/Mjj (which particular mass is used in the denominator does 
not matter in the approximation considered); the normalization V is defined as 

qOmax (PY 

d^ dq °- (150) 

This determines the value of fi 2 normalized at point \x. 

Again, similarly to the situation with A, one can get a lower bound on /i 2 without 
waiting till measurements of the double differential distributions are done. To this 
end one combines the third sum rule ( |149| ) with the Bjorken sum rule and gets [51 



Hi > 3Al ^p 2 - i) (151) 

where the quantities on the right hand side are the same as in Eq. ( |148| ). Numerically 
the right hand side is close to 0.5 GeV 2 , see Ref. |51 . 



6 Quantum-mechanical interpretation 



A quantum-mechanical approach to the derivation of the sum rules for the heavy 
quark transitions has been suggested by Lipkin [27]. The formalism he exploits is 
similar to that used in the theory of the Mossbauer line shape; it is discussed in 
detail in Lipkin's text-book [52| . 

Some of our results can be readily understood within the framework of Lipkin's 
approach, and, therefore, it seems instructive to provide a dictionary allowing one 
to translate (where possible) the field-theoretic consideration in the language of 
quantum mechanics. 

As a matter of fact the expressions for the moments of the distribution functions 
(i.e. averages of powers of the excitation energy) obtained in Sect. 4.2 have very 
transparent quantum-mechanical interpretation, which was already mentioned in 
brief. Namely, the b —>■ c transition in the semileptonic decay is an instantaneous 
replacement of the b quark by c: 



\*c) 



d 3 xe iqx c{0,x)Tb(0,x)\B) = jJB) 



(152) 



where T is some Dirac matrix (7^ for the vector current and 7^75 for the axial 
one). Using the fact that both b and c quarks are heavy we can use non-relativistic 
expansion for the current j^. 

For definiteness, we will consider here the zero recoil limit which was discussed 
above in most detail. Generalization to q 7^ is straightforward. As an example, 
for the vector and axial-vector currents we get 



tf<?= o) = / A (*, [1 - i (i. - i-) 2 + o {£) 



Mx) (153) 
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j£(q = o) = fd 3 x cpl 



IX) 



Ok 



5m; 



[an a k 



imt 



;a k {an) + 



+ 



1 



+ o{ — 



(154) 



where (p c ^ are two-component nonrelativistic spinor fields defined in Eq. (20). Other 
components of the currents do not contain the leading (l/m)° terms. Eqs. ( 153|) , 
( |154j) are still written in the second-quantized form where (p Cj b are field operators. 
However, since the number of heavy quarks does not change in the process at hand 
it is convenient to proceed to the first-quantized form. It corresponds to an ordinary 
nonrelativistic description by a two component wave function 



^afe {flight}) , 



1,2 



where xq is the heavy quark coordinate, a is the heavy quark spinor index and 
{^nght} represents an infinite number of light degrees of freedom. Note that it does 
not imply any non-relativistic approximation for the light cloud; {a^ght} are still 
field-theoretical coordinates of QCD. From now on we will not write out explicitly 
these coordinates in the argument of the wave functions. 
Using these notations Eq. (|152|) takes the form (at q = 0) 



* (Xc W {•••}) 



O 



77T 



*< B >(s , {...}) (155) 
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Jo, and 



(Xc), 
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,m 6 



Jk 



(156) 



- , , , * (B) (x Q , {...}) for j 

Am c mb Km 6 /} 

The derivative appearing in the momentum operators n l = —i{d/dxQ)-\-A l (xQ) acts 
on the coordinate xq of the wave function ty a while a matrices act on its spin index. 

In relativistic and nonrelativistic theories the expectation values (B\...\B) are 
normalized differently. In relativistic theory (B\B) = 2MbV (in the B meson rest 
frame) where V is the volume of the "large box", while (v[/( B )|v]/( B )) = 1; the degree 
of freedom associated with motion of a system (B meson) as a whole is usually 
not considered being treated separately. Correspondingly, the transition rule is as 
follows: 



B(p)\ / d 3 xbOb(x)\B(p)) sccond 



quant. 



(VW\0\¥V) m (157) 



7Tj, etc. 



where O is some local operator, e.g. 70, 

The characteristic feature of the currents considered above is that the states 
produced by them from ty( B ^ are close to \p( D ) and \I/( D *) « a k ^( D \ respectively. 
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Moreover, at q = the currents do not contain terms linear in 1/m; corrections 
start with operators that explicitly contain 1/m 2 . Therefore, 



l- 1 -(---)\B\(anf\B) 



¥ B \x Q ) 
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(o-vrj-o-fc - — -^a k (an) 2 + (an)a k (an) \B)+0 

8mj Am c m b J 

(159) 

Here \B n ) are excitations of B meson, vE^- 8 ™) are the corresponding wave functions. 
(Here and below in this section all matrix elements are written down in the nonrel- 
ativistic normalization; the factor (2Mb) -1 characteristic of the relativistic normal- 
ization is absent.) 

Let us emphasize that v|/( Bn ) are the eigenfunctions of the Hamiltonian 

l-iimo) with ijiq = m b , and not the one with mg = m c . In the final state the 
observable hadrons are described by the eigenstates of 7i(mq = m c ). Therefore 
the amplitude of an individual state \D n ) is obtained by projecting |\I>( Xc )) , \^ k Xc ^) 
onto these eigenfunctions. The projection onto the ground state is equal to unity up 
to 1/m 2 terms, whereas the excited states are produced with the amplitudes which 
scale like 1/m. This reexpansion is the only source of 1/m terms in the amplitudes. 

With this information at hand let us return to the sum rules. We start from the 
first sum rules (see the expressions for Iq V , I^ AA in Eqs. ( |1 1 1| ) , 



l , ± in 1^40. yum , ( |107| )). These sum 
rules have the meaning of the total probability of hadronic production at zero recoil if 
it is assumed that arbitrary energy can be carried away by the lepton pair produced 
in the decay process (see Eqs. ( |116|) , ( |114| )). In the nonrelativistic language these 
total probabilities are nothing but the norms of the states fy( Xc ^ and ^jj. . There is 
no need in reexpansion in order to find these norms. As a matter of fact, the result 
is given by the coefficient in front of vpC 8 ) (or a k ^^ B '); other terms contribute only 
at the level of 1/m 4 . The concrete expressions, thus, are 
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(161) 
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which exactly coincides with the results presented in Eqs. ( |116| ), ( |114| ). 

It is natural to ask at this point how could it happen that the relativistic current 
by^c or 67^750 that was seemingly normalized to unity (through equal time commu- 
tators) lead us to nonrelativistic current with the normalization different from unity. 
The answer is that in the nonrelativistic consideration one excludes the states whose 
excitation energy ~ itlq. In fact, adding such states would restore the normalization 
back to unity. This explains, in particular, the negative sign of 1/m 2 corrections 
to the nonrelativistic normalization. (Let us mention in passing that Lipkin in his 
analysis did not consider these 1/m 2 relativistic corrections altogether.) 

The fact that the states ^^ b \xq) and ak^^ B \xo) are not the eigenstates of the 
Hamiltonian 7i(m c ) becomes important for higher moments I n , Eq. (|106|) , which in 
the quantum-mechanical language take the form 

l Y n v = (¥ B \x c )\ {H(m c ) - (D\H(m c )\D)} n \¥ B \x c )) (l + O (i)) (162) 



and 



1 



In A = \(^ {B) (x c )\a k {H{m c ) - (D*\H(m c )\D*)} n a k \¥ B \x c )) (l + O 

(163) 

The 1/m 2 corrections to currents do not show up in the leading term presented 
above. From now on we will abbreviate (D\Ti,(m c )\D) = (7i(m c ))£>, and likewise for 
D*. 

As we have already seen, all excitation probabilities are proportional to 1/m 2 
and, therefore, 

\n+2 
iv QCD 



In for n > 1 

mi 



The very same 1/m 2 behavior is also readily seen from Eqs. ( |162|) , ( |163| ). Indeed, 
H(m c ) — {7i{m c )) acting on the states (x c ) , a k ^ B ' (x c ) produces a nonvanishing 
result only at the level 1/m: 

(H(m c )-(H(m c )) D )\^ B \x c )) = 

= (H(m b ) - (H(m e )) D + H(m e ) - H{m h )) \^ {B \x c )) = 

((M B - M D ) - K - me) + l* (B) (ay> + O (^) (164) 

and, by the same token, 

(KW - (HWW fft|* (B, (f c )) = K ((M B - M D .) - (m b - m c ) + 



m 2 



Using these equations it is easy to obtain the predictions for the moments 
with n > 2 . To this end we act by the left-most and the right-most operators 
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7i(m c ) — (H(m c )) onto (ty^] and l 1 !/^) , respectively. This produces an overall 
1/m 2 suppression. For vector currents we have 

^ (<*>" (i - i) + *) (« " ^ (<**>" (i - i) + 

+ > ( 166 ) 



m 2 



while for the axial-vector currents 



In A = l^ (B) Wk (H(m c ) - (H(m c )) D *) n a k \¥ B ^ = 
f (W (J- - -LJ) + ^ a fc + -U fe , aB] ) (H - (H))^ 



3 VV \2m r 2m/,/ / 2m, 



x 



x (•> (< w > 2 (i - i) + - ^ ) i* ,B, > + ° (i) (167) 

where 

5 y = (M B - M D ) - (m b - m c ) , 5 A = (M B - M D *) - (m b - m c ) 

are given in Eqs. ( p. 12 ), ( |105| ). Note that in the remaining operators (7i — {7i)) n ~ 2 
one can already neglect the explicit heavy quark mass dependence. For n > 2 the 
terms with 8 lead only to corrections higher than 1/m 2 and can be omitted. Indeed, 
picking up 5 rather than (c?7?) 2 allows one to act by (TC — (H)) directly on \l/ or 
cTfc^ which results, in turn, in an additional 1/m suppression. The case n — 1 is 
somewhat special and will be discussed shortly. 

One can easily see that Eqs. ( |166| ), (|167|) coincide with our previous results for 
the moments, Eqs. (|111|) , (|107|) , provided that one identifies ttq 1 ' 2 ^ in the latter 
with (7i — (H))^ 2 ^ in the quantum-mechanical relations. This correspondence has 
a transparent physical meaning and can be readily proven. Indeed, both play the 
role of the generator of the time evolution, (T~L — (7i)) in the first-quantized approach 
and 7To = ido + Aq in the second-quantized formalism. 

To see that this is indeed the case, let us consider the state 

|*o(t)) = b(x = 0,t)Ob(x = 0,t)\B) (168) 
where O is an arbitrary operator. Then 

(tt o (0)|tt o (*)) = (B\b(x = 0,t = 0)O + Te- i ti A ° { * =0 > T)dT Ob(x = 0,t) \B) (169) 

in the leading in 1/m approximation when the heavy quark Green function reduces 
to the T-exponent above. Expanding the right-hand side in t one gets 



OMO)|*o(*)> = (B\b(x = 0,t=0)O+ J2(ido+Ao) n ^-Ob(x = 0,t)\B) 



(170) 
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that proves the assertion above. 

Let us return now to the discussion of the special case of the n = 1 moments (the 
second sum rule). This case is singled out because, on one hand, Ji must be pro- 
portional to AQ CD /m 2 on general grounds, and, on the other hand, this suppression 
does not immediately show up as it happens for n > 2, since now we have only one 
power of 7i — (H) sandwiched between the bra and ket states. To calculate l\ v and 
I AA we need to account for all 1/m 2 terms that have been neglected in Eqs. ( |164| ), 
( |165| ). Including these terms we get 

lX v = ((M B - M D ) - (m b - m c )) + (¥ B \x c )\ (JL - JL) (an) 2 + 

4 




m 2 in 



^)f-(DE) + 2a-Exn)\¥ B \x c )) + o(^^) , (171) 



and 



I? A = ((M B - M D *) - (m b - m c )) + (^ B \x c )\ - ^-j (na) 2 - ^-BB + 

- -V) (-(DE) + a-Ex n)--l-a-Exn\¥ B \x c ))+0 (^f) (172) 

where we have used the expression ( pi]) for the 1/m 2 part of the Hamiltonian. 

If one neglects the 1 / m 2 terms then l\ v , I AA vanish, and the known relations 
Eqs. (|112|) , ( |105|) for hadron masses at the level 1/m arise. Eqs. ( |171|) , ( |172| ) are 




extension of these relations to the order 1/m 2 . There are two sources of 1/m 2 
corrections in the right-hand side: the first one is explicit 1/m 2 terms in these 
equations; the second source is an implicit (1/m 1 ) dependence of the expectation 
values of operators n 2 and aB over the state ^^(xc). This dependence is due to 
the fact that ^f^ B \x c ) is the eigenfunction of the Hamiltonian T-C(mb) rather than 
of the asymptotic one 7i(mq = oo). It is obvious, by the way, that these matrix 
elements do not produce l/m c dependence. 

Let us discuss in brief how this quantum-mechanical formalism works in the case 
when q ^ 0. The instantaneous b c transition now not only replaces b by c but 
also boosts the c quark providing it with the spatial momentum —q. As a result the 
wave function of the system produced takes the following form for the vector and 
axial-vector currents, respectively: 



¥ x *\x Q , {...}) = e"^ (l - y) ■ {•••}) + 



+ [v\ , (173) 
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+ o (v\ ^£1, ^] (174) 
\ m m 2 J 

where v = —q/m c is the velocity of the final c quark in the SV limit. The factor 
1 — v 2 /8 reflects the overall normalization of the currents; it results in the factor 
1 — v 2 /4 in the n = moments (the Bjorken sum rule). In higher moments it is not 
essential at the level of accuracy we accept here. 

The moments I n {<l) are defined in Eq. (|129|) , but here we define e for the vector 



case as e = M B — q — y + q 2 . The generic formulae for the moments take the 
form: 



I™ = I ^ = 1- V T , (175) 
/ z \ n / A ^ \ 
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jVV = ^(B)| ( e WQ<H( mc ) e -WQ - ^M 2 D + q ^ \¥ B \x c )) + O 

'(176) 

I AA = ^(¥ B ^\a k (e^H(m c )e-^ - y/ M*. + q 2 )" a k \¥ B \x c )) + (177) 



+ o(v\ A2qcb 



v m 2 J 

Let us notice that the only effect of the exponent is to replace the operator 7? in the 
Hamiltonian by 7? — q : 

e*®*H(m e ; n)e~ iq ^ = H{m c ; 7? - g) = 

= ^( mc ; 7 T)-^ + il--l^g. ( TxE + of^ . (178) 

As an example, let us consider the g 2 dependence of the first few moments. For 
n = the result for dlo/dv 2 at {7 2 = coincides with the Bjorken sum rules; it has 
been derived in Sect. 5 with even better accuracy. 

The next moment to consider is n = 1. Here we get for dl\/dv 2 at v 2 = 

dl vv dI AA A 

k= = k=o =2 ' A ~ Mo — m c ~ Mq» — m c . (179) 



This is Voloshin's "optical" sum rule |24 . 

Finally, let us mention n = 2 case, where at v 2 



jrW j tAA 1 

^ ,„ = ^ ,„ = i <0 . (180) 



Concluding this section, let us make a few comments on those aspects of the 
quantum-mechanical approach which we modified compared to the original Lipkin's 
presentation. The central point of the framework suggested in Ref. [37] is the fact 
that the wave function of the charmed system immediately after the instantaneous 
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b — > c transition is known in terms of the B meson wave function, see Eqs. ( |158j ), 
( |159|) and ( |173| ) , (|174|) . These equations differ from their counterparts in Ref . |27j in 
the normalization of the currents; there it was effectively set equal to unity. In our 
expressions the charm wave function is not normalized to unity, with the corrections 
appearing at the level 0(AQ CD /m 2 ) or 0(v 2 ). This correction to normalization 
affects only the first sum rules, viz. Iq V , Jq . 

Another point of difference is that Lipkin did not account for the fact that the 
spin part of the Hamiltonian depends on vtlq at the same level as the kinetic energy. 
As a result, our sum rules at zero recoil are, strictly speaking, different from those 
of Ref. [p7f| . (Our results coincide with those of Lipkin provided the spin interaction 
is switched off.) At the level of accuracy we accept in this section the spin terms 
omitted in Ref. |27J do not affect at all the sum rules at q ^ 0. In particular, the 
linear in q 2 part in the first moments (see Eq. ( |179|) for l\ v and lf A ) derived in 
Ref. 37J identically coincides with Voloshin's sum rule [E3J, although surprisingly it 



was not recognized in Ref. [27]. 

Our final remark concerns cb currents which vanish in the nonrelativistic limit. 
For example, in Sect. 4.4 we considered the spatial components of the vector current 
crfkb at q — 0. In the nonrelativistic limit this current takes the form 
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This current produces the charmed state with the following wave function: 
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The normalization of this state is 
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This expression is identical to the left-hand side of the sum rule ( |121| ) which is the 
n = moment for this current. For higher moments we get 

4 1)VV = li:(*k(xQ)\ (H(m c ) - (H(m c ))) n \M*q)) * 
fc=i 
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\2m c ' 2 mb J* k+ {^r^J t[ ° X M l¥B)) ' (184) 
Now the suppression 1 / m 2 is explicit - it comes from the current itself, and all 1/m 
terms in 7i(m c ) can be omitted. 
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7 Impact of the Perturbative Evolution of Oper- 
ators in OPE 



In this section we shall briefly discuss practical modifications that arise in the cal- 
culation of nonperturbative corrections in the heavy flavor decays if one accounts 
for the normalization point dependence of the corresponding operators. Although 
this question is rather standard, we feel the need to dwell on it in view of apparent 
confusion taking place in applications of the heavy quark expansions existing in the 
literature. 

We have already mentioned above that HQET is nothing else than a version of the 
Wilson operator product expansion. It is well understood (see, for example, a recent 
discussion in Ref. [|28|0 that consistent incorporation of nonperturbative effects in the 
framework of the Wilson procedure requires separating momenta below and above 
a certain normalization point /i. The low-momentum physics is then attributed to 
the matrix elements of the operators whereas the high momentum part enters the 
Wilson expansion coefficients; both, therefore, depend explicitly on /i in such a way 
that all observables are /i independent. 

This procedure is always performed when the corresponding operators undergo 
logarithmic renormalization, for an obvious reason: the Feynman integrals deter- 
mining the coefficient functions in this case logarithmically diverge in the infrared 
domain, and one merely cannot put /i, the infrared regularization, to zero. In 
calculating the power corrections to the heavy flavor decays another situation can 
typically arise - when the operators under consideration have vanishing anomalous 
dimensions and no logarithmic mixing. This is the case, for instance, with the lead- 
ing operator QQ and the kinetic energy operator Qif'Q (the second one still mixes 
with the first one non-logarithmically, through power of /i 2 , see below). In this case 
the coefficient functions possess a safe infrared limit, and it is very tempting just 
to calculate them in perturbation theory per se, with no infrared cut off (or, which 
is the same, putting /x = 0. ) This is what is routinely done with respect to the 
coefficient of the operator QQ. 

From purely theoretical point of view there is no way one can justify such a 
procedure. It gives rise to questions which have no consistent answers, e.g. how the 
sum of all perturbative terms is to be understood, etc. One of the manifestations 
of these inconsistencies is the infrared renormalon (see a recent discussion in Refs. 
|28| , p9|j and references therein). It should be very clearly realized that in the consis- 
tent operator product expansion (and, hence, in HQET) one does not discriminate 
perurbative contributions versus nonperturbative, but, rather, large-distance ones 
versus short distances (all distances are measured in the scale of /i^ 1 ). 

However, in practice no separation of the infrared part from the coefficient of 
QQ is carried out. The usual routine is as follows: one takes the known expressions 
for the one-loop (or two-loop) perturbative corrections for a particular quantity and 
merely adds to these perturbative terms nonperturbative corrections expressed via 
certain matrix elements. (We also follow this routine for numerical calculations). 
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For instance, the perturbative correction r/A to the the axial-vector current 07^756 
calculated from the standard Feynman graphs with no subtractions was simply 
added to the nonperturbative contribution of Eq. ( |115[ ). It is clear however that 
the nonperturbative contribution per se takes care of all relevant gluon exchanges 
with momenta below fi; on top of it the one-loop Feynman integral for the radia- 
tive correction t]a has some (small) piece coming from the same low-momentum 
domain. The question which immediately comes to one's mind is the menace of 
double counting. 

The answer to this question is that in actual OPE-based calculations, which are 
always truncated at some finite order in a s and keep only a few higher-dimensional 
operators, one can stick to what is called the "practical version" of OPE in QCD 
and still avoid double counting. 

Indeed, it is not possible to define the perturbative part in, say, (Qt? 2 Q) to 
all orders in a s . On the other hand, it is quite possible to introduce a "one- loop 
perturbative contribution" to Qn 2 Q normalized at fi. To this end we, by definition, 
take two gluon lines in Qn 2 Q, contract them to form the gluon loop, use the bare 
gluon propagator and calculate the loop cutting the integration off from above, at 
kg = \i. We then get 

2 2 

Qll Q\onc— loop T, A' QQ ■ 

Now let us subtract and add this "one-loop" 7r 2 from the properly and scientifically 
defined Qn 2 Q. Moreover, let us combine the added part with C(fi)QQ; then we get 
QQ times the coefficient coinciding - up to corrections of higher order in fi 2 /m 2 - 
with the coefficient obtained from the full perturbative one-loop calculation, with no 
subtractions whatsoever. Simultaneously, the matrix element of the kinetic energy 
operator is replaced by that with the subtracted "one- loop" part. Strictly speaking, 
this quantity does not represent now matrix element of any operator; hence, fac- 
torization inherent to the genuine OPE is lost. This is unimportant for numerical 
analysis due to the fact that the term which we added and subtracted is numerically 
small, much smaller than the actual value of (vr 2 ). For this reason the replacement 
of the added term plus C{n)QQ by merely C onc -i oop QQ introduces a very small 
numerical error, much smaller than it 2 correction itself. This fact - the numerical 
smallness of the "one-loop" value of the condensate compared to its actual value - 
constitutes the conceptual foundation for the practical version of OPE. In this ver- 
sion fi does not appear explicitly in the coefficient function of the lowest-dimension 
operator; formally this corresponds to setting fi = there. The fact of the numeri- 
cal smallness is not an obvious property of QCD and is to be cross-checked in any 
new situation. So far it turns out that it is always valid, for reasons which are not 
completely understood, see the reprint volume cited in [[53|. We repeat, however, 
that in principle one should calculate the Wilson coefficients by evaluating relevant 
Feynman integrals with an explicit infrared cutoff in the propagators. 

To further facilitate understanding of this subtle aspect let us discuss this general 
strategy, which we usually stick to, in a particular example, namely the first sum 
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rule, Eq. ( |114| ). For simplicity we will forget about the operators that depend on 
spin of the heavy quark since they do not mix with the leading operator and are 
irrelevant for the aspect under discussion. The sum rule ( |114j ) with the perturbative 
corrections added takes the form 

Fb-,D* + ^B-»excit = U ~~ oZ^2 ~ ^I—J^L ( _ _| J ! (185) 

i= Y^ 3m~ 4 \m l c m{ 3m c mb J 

nonperturbative corrections are then represented only by the kinetic operator 7F 2 . 
For a consistent calculation of the perturbative correction factor U we need to 
introduce an infrared renormalization point fi -C m c and use Ufa) m Eq. ( |185| ); the 
operators entering this sum rule are then normalized at the scale fi. It is worth 
noting that without introducing fi the sum rule written above, strictly speaking, 
has little sense because the sum over excitation will diverge in the ultraviolet (see 
below). If a renormalization point is set, the sum will run only over states with 
excitation energy below fi. Keeping in mind the explicit /i-dependence of Ui it is 
clear that it cannot be equal to i]\ which is calculated without an infrared cutoff 
and is thus /^-independent. In reality, however, the major part of both rf A and Ufa) 
comes in the heavy quark limit from the momenta ~ rriQ and therefore they must 
be similar. Actually, one even do not need to calculate the perturbative factor U 
anew as long as r\\ is known. 

To determine U one notes that to any particular order in perturbation theory 
the relation holds 

r)\ = Ufa) Im=o ! (186) 

in fact it is merely the definition of t\a- It can be formally obtained considering the 
sum rule (|185|) in the perturbation theory in the particular order (we will assume the 
first order in a s in what follows). In this approximation the spectrum of particles is 
given by quarks and gluons, and matrix elements of operators, in particular are 
nothing but "perturbative one loop" contributions considered above. The elastic 
contribution is given by the exclusive b — > c probability whereas the sum over 
excitation is represented by the process b — > c + g whose probability is proportional 
to a s . Then, at fi = one gets the relation ( 186 ) because other terms vanish. 



To determine /1 dependence of U we can use the same technique as in Sect. 4.4, 
namely, consider the perturbative analog of the sum rule (|185|) : 

+ ^ry 3 j^—\(cg\c^ 5 b\b)\ 2 = 
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The gluon emission amplitude is now 
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Here n = k/\k\. Using the expression ( |126| ) for (/i 2 ) pert we immediately find that 

U ^ + ^_ + ^ + _] . (189) 

Similar situation occurs in the sum rule for the vector current (the time component) 
where both the coefficient in front of /i 2 and the gluon emission probability are 
proportional to (l/2ra;, — l/2m c ) 2 . 

The scale dependence of the strong coupling a s in the above expressions is not 
essential in the one loop calculations; if higher loops are accounted for, the sum of 
all Feynman integrals with the particular infrared cutoff fi will automatically have a 
suitable form to give logs of the ratio of mj \x necessary to convert a s {m) appearing in 
the perturbative calculations of £, into a s (/x) in terms responsible for /i 2 corrections. 

In principle, one can now utilize Eq. (|189|) to use the exact perturbative value of 
£a(aO in Eq- ( |185| ), which differs from rj A for /i ^ 0. Obviously, ^(//)— T] A in Eq. (|189|) 



contains two pieces, one given by the perturbative one loop sum over excitation 
calculated with the upper cutoff /i, and the second representing the "perturbative 
one loop" contribution to fJ^(fj)- It is clear that one can formally subtract these two 
terms from the sum over excitations in the left hand side of the sum rule ( [L85|) and 
from /i 2 (/i) term in the right hand side, respectively, and then use as the perturbative 
factor £a hs one loop value at fi = 0, i.e. r) A . It is just what one does following the 
routine practice of OPE where the perturbative coefficient is calculated without an 
explicit infrared cutoff. In such a case /x 2 entering the "practical" form of the sum 
rule can be merely understood as 



and, simultaneously, 



A^^M-^A (igo) 



</i 2uj(2tt) 3 Am c mh 



The normalization point fi is chosen high enough to ensure the approximate duality 
of the parton computed probabilities to real hadronic ones above \i. Then one can 
remove the upper limit in the difference in the right hand side of the above equation 
and assume that the sum runs over all energies: the region above \x cancels out 
automatically. 

The expression in the right hand side of Eq. ( |190| ) is nothing but a linear (in 
/i 2 ) extrapolation of the matrix element from the point \x to fi = 0. Formally, it is 
independent on /i : the dependence appears in terms proportional to a 2 s . Therefore 
one can take any value of \x as long as Aqcd *C fi m c . At the same time, the sum 
over excitations now is to be understood, strictly speaking, as the one from which 
the small one loop perturbative gluon probability is subtracted. 
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A similar analysis can be literally extended to the two, three, and in principle, 
to any finite loop calculation of £4. We would like to stress once again, though, 
that this procedure has no theoretical justification and does not carry any physical 
meaning. We dwell on it only for the reason that it literally corresponds to the 
routine procedure used in numerical calculations; using /x = not only does not 
allow to use consistently OPE, but would lead also to numerical problems in higher 
orders. 

Now let us turn from this rather general theoretical discussion to more practical 
questions related to heavy flavor decays. Up to now nonperturbative effects have 
been discussed in detail through corrections of order I/tuq. Some effects, like the 
invariant mass of the final hadronic state in the decays considered above, have 
corrections starting at order 1/ttt.q and are expressed via the parameter A (see also 
Ref. |fLl|| ). As pointed out in Ref. |28| and illustrated in the present paper, their 
effects are determined by Feynman integrals which have a linear behavior in the 
infrared region; the corresponding IR effects are not expressed in terms of matrix 
elements of any local operator. 

The inclusive widths of heavy flavor particles are obtained from an expansion 
in local operators, and corrections start with terms scaling like 1/wJq. These are 
described by two universal operators - the chromomagnetic one, Q ^aG Q, and the 
kinetic operator, Q (iD) 2 Q. The natural normalization scale for them is given by 
fi ~ itlq, at least if one neglects the mass of the final state quark(s) One cannot 
however use directly this high normalization point because then the matrix elements 
of high dimension operators will scale like itlq to the corresponding power due to 
purely perturbative contributions, and instead of an expansion in 1/mg one would 
obtain the suppression of the higher order terms only as some powers of a s (mo). 
To obtain the real power expansion one must evolve these operators down to a scale 
/i which is to be much smaller than tuq but still much larger than Aqcd- The 
expansion one arrives at in this way runs, strictly speaking, in powers of h/itiq. 

The present state of the art in this kind of calculations is limited only by one 
loop corrections to Wilson coefficients which - apart from the chromomagnetic op- 
erator that has a logarithmic renormalization - are calculated (or even typically 
borrowed from old QED calculations) without an infrared cutoff. The same refers 
to the corrections to weak currents used in the present paper. The analysis above 
suggests therefore that the value of the kinetic term /i^ can be understood in the 
corresponding expressions as a linear extrapolation to fi = 0. This problem does 
not arise at all for the chromomagnetic operator Oq whose mixing with the leading 
one, QQ, appears only in the next order in I/tuq due to the fact that it is not a 
spin singlet; no double counting occurs for this reason. From a practical viewpoint, 
because the value of /ijj. is basically unknown yet, it does not make a big difference 
at present to prefer this or an alternative definition. Let us note in passing that it 



15 It is important that for the hyperfine splitting in heavy mesons which allows one to extract 
experimentally the chromomagnetic matrix elements, the same normalization point emerges. 
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is quite probable that the QCD sum rule estimates determine a similar quantity ex- 
trapolated to n = because no explicit infrared cutoff in the integrals is introduced; 
though this question definitely deserves a more careful analysis if one wants to go 
beyond the accuracy of "practical" OPE. 

The issue of an accurate understanding of the definition of matrix elements of 
operators becomes important when one turns to the real practical bounds on physical 
observables of the type discussed in this paper. For example, extrapolating the 
operator Q tt 2 Q to the zero renormalization point implies a subtraction of a positive 
quantity which, in principle, might have even changed the sign of the matrix element. 
To state it differently: one may be concerned whether the inequality (|1 23|) survives 
the extrapolation of fi 2 to a low point, which is often assumed implicitly. We shall 
argue now that this effect is too small numerically and cannot upset the bound we 
used. 

To see it, let us consider the reasonably high normalization point /j, ~ 1 GeV. Us- 
ing the explicit estimate of the renormalization point dependence of operator Qn 2 Q 
in Eq. ([7^) and assuming a s (/z) ~ 0.36 one readily obtains that the amount one may 
need to subtract from fi 2 (fi) constitutes at most 0.15 GeV 2 , a value that does not 
exceed the theoretical uncertainties in the existing estimates of fi 2 . Most probably 
this number overestimates the real contribution to be subtracted, because approxi- 
mate duality of the perturbative corrections is expected to start earlier; moreover, 
the perturbative coefficient functions (tja, rjvi corrections to inclusive widths etc.) 
are evaluated in one loop using a smaller value of a s ~ a s (mg). The second effect, 
though formally of higher order in a s , is too transparent physically to raise doubts 
that a more realistic estimate of effects of potential double counting corresponds to 
using a s {m c ,m h ) rather than a s {l GeV) above. 

At the same time, as emphasized in Sect. 4.4, if the normalization point is intro- 
duced via the upper bound in the integral over the energy of the excited states, the 
inequality fi 2 > fiQ holds for any normalization point (in other words, such regular- 
ization does not violate the positivity of the Pauli operator for the spinor quark). 
Therefore, it is legitimate to take the normalization point as low as 1 GeV. In this 
case, obviously, one deals with fj,fj normalized at this low point as well, and it is 
known ]54|| that the perturbative evolution increases its value toward lower ! In 



our estimates we took /i^ directly from the hyperfme splitting of B and B*, therefore 
that value corresponded to /i ~ 4.5 GeV. Apparently its hybrid log enhancement 
would safely make up for the relatively insignificant subtraction of the "perturba- 
tive" contribution to fi 2 . Based on these arguments we have stated in the previous 
paper HI| that the inequality 



2 ^ 2 



must survive in QCD, in spite of recent claims |49| that it cannot hold true. 

It is instructive to trace how this inequality works at different scales [i. Most 
trivially it is fulfilled when /i is taken parametrically large. Then fi 2 contains large 
positive perturbative piece of the order of — fi 2 that grows faster than any possible 
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change in /Iq having no additive renormalization, even if the hybrid anomalous 
dimension of the latter were negative. 

A more interesting consideration emerges when one wants to push /i toward lower 
values. Using the naive one loop expression for the evolution of /j/q corresponding to 
the hybrid anomalous dimension 7g = 3 one would obtain an arbitrarily large value 
for which, of course, makes little sense. The answer to this apparent paradox 
is rather obvious, especially if one looks at a hypothetical zero recoil excitation 
curve (for the external current cjib ) similar to the one depicted in Fig. 4. The real 
evolution of the difference fi^ — ji 2 G at m q 3> uiq, according to the sum rule ( |121| ), 
is given by the decay probability occurring at energy e = /i; obviously the latter in 
no way is given at low e by the simple perturbative formulae based on the strong 
coupling a s with the Landau pole, and rather stays finite at any ji. In other words, 
the evolution of fiQ is to be smooth even when one approaches the strong coupling 
regime, and no formal contradiction emerges. 



8 Conclusions 

In the present paper we have addressed weak transitions between heavy quarks from 
the "inclusive" side most suitable theoretically for applying the technique of the 



Wilson OPE. This analysis naturally extend the consideration outlined in Ref. |I5 
which concerned the heavy quark distribution function relevant for the decays in the 
limit of small velocity for the final state hadron system. It has been demonstrated 
that a few sum rules discussed so far in the literature are in fact dispersion relations 
for the moments of a single distribution function, considered in different orders in 
1/tuq and in different kinematical regimes. We have shown how the expansions 
of HQET can consistently be obtained from QCD using this strategy, and in this 
way illustrated that such natural assumptions as "global duality", that usually are 
attributed only to the inclusive width calculations, are in fact necessary ingredients 
in any consistent model-independent treatment and, in particular, are implicitly 
used in HQET as well. 

The analysis of the sum rules proved to be very instructive in elucidating the 
important fact that has usually been neglected in HQET - the necessity of introduc- 
ing an explicit infrared normalization point n ensuring true separation of low and 
high momentum physics, which cannot be set to zero. The consistent application 
of this approach leads to the fact that all nonperturbative parameters, including 
A, cannot be sensibly defined as universal constants, but rather depend explicitly 
on the normalization point. This has been previously mentioned in our paper [15| 



and discussed in detail in Refs. |28], |29j. Here we gave a physical illustration of 



how it works analyzing possible constructive phenomenological definitions of corre- 
sponding quantities in the presence of radiative corrections. In this way we have 
supplemented the previous calculations by estimates of the dependence of the kinetic 
energy operator Qif 2 Q on the renormalization point. 
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The fact that such "purely nonperturbative" objects like the pole mass of the 
heavy quark, A, "purely nonperturbative" distribution function of heavy quarks 
routinely used in HQET are incompatible with a consistent OPE-based approach 
and are ill-defined theoretically calls for the clarification of how the known results 
on nonperturbative corrections in HQET must be interpreted. This does not mean 
of course that the concrete calculations that have been done so far are irrelevant, 
and we formulated the way in which they are to be understood for a few typical 
examples. 

As a practical application of our sum rules we have derived a model independent 
lower bound on the deviation of the exclusive axial form factor Fb^d* of the B — > 
D* + I v decay at zero recoil - a process that for a long time has been believed to 
give the best theoretical accuracy to determine and estimated a reasonable 

'central' value, F B ^ D *(q = 0) ~ 0.9 . The deviation appears to be essentially larger 
than the estimates that had been obtained before from model calculations based on 
standard HQET, and apparently better agree with quite general expectations about 
the size of corrections to the Heavy Quark symmetry for charmed particles. On 
the other hand, the theoretical clarification of the notion of the heavy quark mass 
made in recent papers J28|, ^9f suggests that the most accurate theoretically way to 
determine the CKM matrix elements for heavy quark decays is using the inclusive 
semileptonic widths. These results have been reported in Ref. ||T . 



It is worth clarifying in this respect that in our estimates of the exclusive form 
factors of the b — > c transitions we consistently took into account terms through 
order I/tuq and discarded effects that scale like l/m%. The parameter l/m c is 
actually not very small and even the second order corrections are as large as 10% 
here; therefore one can expect sizable relative corrections for real form factors due 
to higher order terms. In particular, this applies to the model-independent upper 
bound for Fb-,d*- Our result is strict in the sense that it holds for corrections 
through terms of order l/m,Q that have been addressed in the literature so far. 

One of the sum rules at zero recoil enabled us to rederive a model independent 
lower bound on the value of the kinetic energy operator in B mesons, 

A = ^ B \™ h \ B )*\ (mi, -mi) 

in a way that clearly showed its physical relevance even in real QCD, and not only in 
the approximate framework of quantum mechanical consideration, as it is sometimes 
stated. 



ACKNOWLEDGMENTS: N.U. gratefully acknowledges the creative atmo- 
sphere and hospitality at Theoretical Physics Institute of University of Minnesota, 
and the exchange of ideas related to the subject of this paper with V. Braun. 
We would like to thank B. Blok, A. Falk, T. Mannel, M. Savage and especially 
L. Koyrakh and M. Voloshin for useful discussions. This work was supported in 
part by the National Science Foundation under the grant number PHY 92-13313 
and by DOE under the grant number DE-FG02-94ER40823. 



65 



9 Appendix 

Let us derive the expression 



A = _ l |(ff|7T»| 2 

2 n 3 E n — Mb 



in ordinary quantum mechanics. We will use the nonrelativistic normalization of 
states below. 

It is easy to see that in simple potential description of the heavy hadron as a 
two body system with the reduced mass m r ~ m sp the sum on the right hand side is 
given just by the half of this mass (which in this case is identified with A). Indeed, 

_ \(B\ifv\n)\ 2 



E n — M 



B 



represents the order v 2 correction to the ground state energy, produced by a per- 
turbation 

5H = 7TV . 



On the other hand 



H + 5H = H(n + v)- 



2 

The eigenvalues of the Hamiltonian given by the first term on the right hand side 
are the same as for the unperturbed one, which leads to the relation 

\- I (B\nv\n) | 2 _ m r v 2 

\ E n -M B ~— {AA) 



which coincides with Eq. (|94]) if A = m sp . In reality, of course, Eq. ( p4\ ) is more 
general and accounts also for the binding energy. To show it we can apply the 
following general consideration. 

First, let us note that the relation similar to Eq. (|A.1|) can be obtained in a more 
general way. Namely, for any system whose Hamiltonian depends on the heavy 
quark momentum ir in the non-relativistic quadratic way, 

7T 2 

H = h H]i S u(xQ, {flight}) (A.2) 

one has the exact commutation relation 

7? = im Q [H,x Q } (A.3) 

where 7? and xq are the operators of the heavy quark momentum and coordinate. 
Then one writes 

_ |(0Hn)| 2 ,mQ / (0|7r fc |n)(n|[^,xg]|0) | (0|[^,xQ]|n)(n|7r fc |0) \ = 
n E n — E 2 n \ E n — E E n — Eq J 
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= ^E ((0|7r fe |n)(n|^i0)-(0|x?|n)(n|7r fe |0) 

* n 

= ^£<Ofe,*?]|0> = ^ . (A.4) 

n 

Note that we do not sum over spatial index k in this equation; it is valid for arbitrary 
k. Below in this section we will assume that the summation over k is not performed. 
We emphasize that this relation is rigorous for any system as long as the heavy 
quark momentum 7? enters the Hamiltonian quadratically. 
Now we apply Eq. ( |A.4j ) to B meson. We get 

m Q _ v-* \( B Wk\n)\ 2 

T~\ E n -M B ■ (A - 5) 

The sum over intermediate states runs over all possible hadronic states which are 
marked, in particular, by their total momentum tt. The matrix elements, however, 
do not vanish only for zero momentum transfer: 

(B(p = 0)\n k \n(p}) = (BW\n) QM {2 * f Jl {P) (A.6) 



where V is volume. When squared only the equal momentum matrix elements are 
present, and the factor (27r) 3 <5 3 (0) = V cancels against l/y/V normalization of states 
associated with the continuous spectrum of total momentum. The situation requires 
more care when the state n is B meson; the energy denominator E h — Mb = p 2 /2Mb 
has a pole in this case when integrated over d 3 p , and this singularity must be treated 
properly. We write therefore 

m Q \(B\n k \n) QM \ 2 f d 3 p \(B(p = 0)\rt k \B(p))\ 2 (A?) 

2 £ B E n — Mb J (2tt) 3 p 2 /2M B * 1 " ' 

The matrix elements in the last term can be represented in the form 

(B(p = Q)\* k \Bffl = ^(B(P = 0)\P k \B(p}} (A.8) 

where P is the total momentum operator. To calculate the integral over the mo- 
mentum of B meson we can again use Eq. (|A.4j ): 

#p \{B(p = 0)\P k \B(p))\ 2 _M B 

(2tt) 3 p 2 /2M B 2 ' 1 ' ' 

this is nothing but the relation (|A.4j ) for "quantum mechanics" of free B mesons 
considered as elementary particles. Combining Eqs. ( |A.7| )-( [AI?D we finally get 

m Q \{B\n k \n) QM \ 2 
2 ^ M n -M B 2M B ' 

or 

, \(B\TT k \n) QM \ 2 = m Q M B - m Q A 
^ M n -M B M B 2 2 1 ' ' 



This equation is clearly the relation ( pm we are looking for; if one sums over k 
additional factor 3 emerges to match the exact coefficient in the latter. 
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Figure Captions 

Fig. 1. The tree graph for the transition operator. 

Fig. 2. A qualitative picture of the spectrum dY/dE^ in the Hq — > (f)X q with 
0(v 2 ) terms included. The monochromatic line of the quark transition Q — > 0g 
(the dashed line at £ = £0) is dual to the physical line corresponding to the elastic 
decay Hq — > H q (p at E = EQ hys plus a shoulder due to the transitions to the excited 
states Hq — > </>i?*. The height of the shoulder is ~ v 2 . Hard gluons are neglected. 

Fig. 3. The diagram responsible for the one-gluon correction in the energy 
distribution 

dT(Q — ► <pq + gluon) 

Fig. 4. A sketch of the energy spectrum dY/dE^ with C(a s ) radiative tail 
included. 

Fig. 5. The cuts of hi(q ) in the complex q plane at q = 0. The point go — 
Mb + Md marks the beginning of three cuts: two of them originate at the mass of 
B c meson and the third starts at M B + M D . The three cuts can be differentiated 
considering q 7^ 0. Only a part of one cut at < q < M B — M D is relevant to the 
decay process under consideration. The masses are not shown in a real scale. 
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